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Abstract 

This  papa  is  the  third  in  a  series  in  which  we  study  the  superconvergence  of 
finite  element  solutions  by  a  computer-based  approach.  In  [1]  we  studied  classical 
superconvergence  and  in  [2]  we  introduced  the  new  concept  of  17%-superconvergence  and 
showed  that  it  can  be  employed  to  determine  regions  of  least-error  for  the  derivatives 
of  the  finite  element  solution  in  the  interior  of  any  grid  of  triangular  elements.  Here 
we  use  the  same  ideas  to  study  the  superconvergence  of  the  derivatives  of  the  finite 
element  solution  in  the  interior  of  complex  grids  of  quadrilaterals  of  the  type  used  in 
practical  computations. 
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1  Introduction 


Classical  superconvergence  studies  the  exceptional  rates  of  convergence  of  values 
of  quantities  (derivatives  of  the  displacement,  strain,  stress  etc.)  computed  from  the 
approximate  solution  at  certain  special  points  in  the  mesh.  These  points  are  called 
superconvtrgtnct  points  for  the  quantity  and  the  corresponding  values  are  called  su- 
perconvergent.  We  distinguish  two  tjrpes  of  superconvergence: 

a)  Direct  superconvergence:  The  superconvergent  values  for  the  quantities  are  ob¬ 
tained  by  direct  evaluation  from  the  approximate  solution  at  the  superconver¬ 
gence  points. 

b)  Superconvergence  via  post-processing:  The  superconvergent  values  are  obtained 
by  some  post-processing  (local  averiq^g)  from  the  approximate  solution. 

In  this  paper  we  will  study  only  direct  superconvergence.  A  similar  study  of  super¬ 
convergence  via  post-processing  is  possible  and  will  be  the  subject  of  a  forthcoming 
paper. 

Classical  superconvergence  for  meshes  of  quadrilateral  elements  has  been  studied 
by  several  investigators;  see  for  example  [3-9]  and  the  references  therein.  It  should 
be  noted  that  all  the  mathematical  studies  of  superconvergence  for  quadrilaterals  el¬ 
ements  have  addressed  the  superconvergence  in  meshes  of  squares  or  meshes  which 
can  be  mapped  to  a  mesh  of  squares  via  a  suiSciently  smooth  transformation  (see 
for  example  [4-5]).  However,  practical  finite  element  meshes  have  to  conform  to  com¬ 
plex  geometries  (boundaries,  material-interfaces)  and  have  to  be  locally  refined  near 
reentrant  comers,  fillets  or  other  critical  re^ons.  In  [1]  we  have  demonstrated  that  the 
superconvergence  points  are  very  sensitive  to  the  geometry  of  the  mesh,  the  solution-type 
and  the  coefficients  of  the  differential  operator.  Hence  the  superconvergence  points, 
in  the  classical  sense,  may  not  exist  for  the  meshes  used  in  engineering  computations. 
Thus  classical  studies  of  superconvergence  do  not  result  in  concrete  guidelines  on  how 
to  sample  the  various  solution  quantities  in  the  practical  grids. 

In  [1]  we  introduced  the  following  new  definition  of  superconvergence:  Let  {u^}  be 
a  one  parameter  sequence  of  finite  element  solutions  of  a  problem  which  are  computed 
using  a.  sequence  of  meshes  T  =  "  denote  the  exact  solution.  Let  us 

assume  that  we  are  interested  in  the  values  of  the  solution  or  its  derivatives  or  linear 
combinations  of  these  quantities  i.e.  in  the  linear  functional  F(u)(e).  Let  us  assume 
that  for  every  element  r  of  the  mesh  a  special  point  x,  which  depends  on  the 
geometry  (but  not  the  size)  of  the  element,  is  given.  Then  denoting 


(1.1)  $(u-uJ:=max|F(u-uJ(a!)| 

we  are  interested  in  the  values  of  relative  error  in  F('u)  at  x, 
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(1.2)  := 


If  the  point  m  is  such  that 


f  l^(«  -  «fc)(«)l 
I  -  «a) 


if»(tt-ttO^0 
if  ^(ti  -  tt,^)  =  0 


(1.3)  0(»;F;ti,i*fc,h,T)  <  ^  as  h — >0 

then  m  will  be  called  a  u-ri%-superconveTgence  point  relative  to  the  exact  solution  u 
and  the  family  of  meshes  T.  Consider  now  a  family  of  solutions  li;  the  point  m  will  be 
a  superconvergence  point  if  it  is  u-i7%>snperconTergence  point  for  every  u  €  U. 

Obviously  @(0;i^;u,u^,h,r)  <  1  and  thus  all  points  in  every  element  r  are  100%- 
snperconvergence  points.  Note  that  if  there  mdsts  a  point  m  in  the  element  r  for  which 
17  =  0  in  (1.3),  i.e. 


(1.4)  0(*;  F;ii,tt;^,h,r) — >0  as  h — >0 

for  a  particular  u  (resp.  for  every  u  E  U)  then  e  is  a  it-superconvergence  (resp.  U- 
superconvergence)  point  tn  the  classical  sense. 

So  far  we  have  assumed  that  we  are  interested  only  in  one  scalar  functional  F.  We 
could  also  be  interested  in  the  points  which  are  superconvergent  simultaneously  for  two 
(or  more)  functionals  e.g.  the  derivatives  in  the  and  x,  directions.  We  can  formalize 
this  by  assuming  a  vector  functional;  the  meaning  of  the  17%-superconvergence  point 
for  the  vector  functional  is  obvious.  Sometimes  we  will  call  it  explicitly  as  simultaneous 
superconvergence. 

Let  us  note  that  if  17  is  very  small  (smaller  than  a  tolerance)  we  can  still  practically 
speak  about  a  superconvergence  point  (with  a  given  tolerance),  which  is  the  reason  for 
introducing  the  notion  of  the  17%-superconvergence.  For  example  in  Section  5.1.1b  we 
report  the  simultaneous  superconvergence  for  the  x.  and  x,  derivative  with  tolerance 
less  than  0.001%  (for  a  mesh  of  biquadratic  elements  with  a  mesh-interface). 

Remark  1.1.  In  this  paper  we  are  interested  in  the  pointwise  superconvergence  as 
defined  above.  Of  course  the  notion  of  superconvergence  can  also  be  understood  in 
different  ways.  For  example  we  can  be  interested  in  the  quantity 

(  E  E  1^(«  -  «/.)(»r)n  * 

(1.5)  #,({®J;F,n,n;„h,{r})  ;=  -  - 

r6{T} 
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where  {r}  is  a  set  of  elements,  {e^}  ^  points  in  each  element.  Then  the 

pointwise  superconveigence  is  related  to  q  =  oo.  Very  often  for  locally  refined  meshes 
the  area-measure  of  nonregular  mesh-structure  is  of  order  0{h).  K  in  addition,  as  it 
often  occurs  (see  Section  5.4),  tf  — »  0  (at  the  points  {0^})  exponentially  with  the 
distance  of  the  element  r  to  the  interface  we  get  — ►  0  as  h  — ►  0  (in  fact  s»  hi). 
Various  authors  are  calling  this  effect  also  as  superconvergence  (see  e.g.  [4,  5]).  Another 
similar  example  is  to  create  an  ‘‘improved”  solution  (see  [10,  11])  and  measure  the 
error  of  this  “improved”  solution  (see  [12]).  Once  more  the  notion  of  superconvergence 
depends  on  its  measure,  e.g.  if  measured  in  Z*-norm,  9  <  00,  we  get  superconvergence 
effects  where  for  9  =  00  there  is  no  superconvergence  (see  the  example  of  a  grid 
with  mesh-interface  in  Section  5.2  below).  Hence  the  word  auperconvergence  has  to  be 
properly  defined. 

The  goals  of  the  paper  are  the  following: 

1.  To  show  that  the  17%-superconvergence  points  in  an  element  depend  on  the  local 
neighborhood  topology  of  the  mesh. 

2.  To  address  in  detail  the  superconvergence  for  concrete  topologies  in  meshes  of 
square  elements  with  refinements.  These  results  can  be  practically  used  for  the 
meshes  which  have  locally  this  character.  (The  meshes  used  in  various  codes 
have  only  a  finite  number  of  local  topologies.) 

3.  To  show  that  if  we  ask  only  for  improved  accuracies  in  an  element  there  are 
areas  where  they  can  be  achieved.  These  areas  are  robust  with  respect  to  the 
mesh-topologies  outside  the  immediate  neighborhood  of  the  element. 

4.  To  demonstrate  that  the  increased  accuracy  can  be  achieved  only  if  the  meshes 
are  such  that  the  effect  of  the  pollution  is  under  the  control. 

5.  To  develop  a  general  procedure  to  address  the  superconvergence  problem  for  grids 
with  general  topologies. 

Following  this  Introduction  we  outbne  the  model  elliptic  problem  (Poisson’s  equa¬ 
tion)  and  we  describe  the  classes  of  meshes  of  quadrilaterals  for  which  the  study  of 
superconvergence  is  given.  We  then  describe  a  computer-based  approach  which  can 
be  employed  to  find  the  points  of  least-error  for  the  derivatives  of  the  solution  in  any 
interior  element  of  any  mesh.  We  give  examples  of  application  of  this  approach  to 
study  classical  and  f/^-superconvergence  for  complex  meshes  of  quadrilaterals  with 
local  refinements. 


5 


2  Preliminaries 

2.1  The  model  problem 

In  this  paper  we  will  study  the  superconvergence  of  the  derivatives  of  finite  element 
solutions  of  Poisson’s  equation.  Let  (1  C  R’  be  a  bounded  polygonal  domain  and  let 
its  boundary  dfl  be  split  into  two  parts  Td  and  Tff  (Fn  has  positive  length).  Let  u  be 
the  solution  of  the  problem 


(2.1a) 


(2.16) 


tt  =  0 ,  on  Fp  , 


(2.1c) 


du  ^  du  _ 


Here  nk,  ib  =  1,  2  are  the  components  of  the  unit  outer  normal  of  dfl;  /  €  L’(n), 
g  €  X^(F;^)  are  jpven  data.  If  F^  =  0  it  is  assumed  that  the  compatibility  condition 

/  o  =  0  is  satisfied. 

Jafi 

Let  Rf^(fl)  :=  |u  €  H^(il)  |  u  =  0  on  Fp|.  Then  the  above  problem  may  be  put 
in  the  variational  form:  Find  u  €  Xrp^(n)  such  that 


(2.2)  Bn{u,v)  =  L„{v)  y  v  e 

where 


(2.3)  Ba(u,v)  ;=  f  Vtt  •  Vv ,  L^iv)  '=  f  fv  +  f  ^  ■ 

JO  Ju 

Let  T  :=  {Th}  be  a  regular-family  of  meshes  of  quadrilaterals  with  straight  edges. 
(It  is  assumed  that  for  any  quadrilaterals  t,-  €  3\,  the  intersection  Ti  D  Tj  is  either 
empty,  a  vertex  or  a  common  edge,  and  that  the  regularity  conditions  (37.40)  in  [13] 
hold.  For  the  regularity- conditions  for  grids  of  quadrilaterals  with  refinements  see  [14].) 
The  meshes  Th  are  not  assumed  to  be  quasiuniform.  We  introduce  the  conforming 
finite-element  spaces 
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(2.4a) 


^:={»€C»(n)l  <€SS(t.), 


*  =  1, 


(2.4t)  Si(T,)  =  {»,  €  C“(r.)|  v,cF^e  5l">(f)} 

where  f  :=  (—1,1)^  is  the  master-element  and  is  the  bilinear  mapping  of  f  onto 
Tj^;  M{Th)  is  the  number  of  elements  in  and 

(2.5)  :=  {  P  I  =  E 

By  our  assumption  the  mapping  is  a  bilinear  map  of  f  onto  Hence 


*1  —  X!  ®»i  ®2  »  ®a  ~  Xrf 

»j=o  »J=0 


Let  P(xi,*j)  =  53  Cron  be  a  polynomial  of  degree  p;  then 

0<m+n<p 


(2.7)  P(li ,  X,)  =  Y,  Cron  (  ^  x\  X^)*"  (  ^  xj  x',)"  =  X!  »'  »J 


0<m+n<p  <^'=0 


Hence  any  polynomial  of  degree  p  belongs  to  the  span  of  the  shape-functions  )  on 
every  physical  element  t,^. 

Below  we  will  consider  meshes  which  are  obtained  from  a  h-refinement  procedure 
which  include  elements  with  irregular  connections.  We  say  that  an  element  has  n- 
irregular  connection  on  one  of  its  sides  if  it  is  connected  to  (n-f-  l)-neighboring  elements 
on  that  side.  Here  we  consider  only  elements  with  1-irregular  connections  and  we  will 
refer  to  the  meshes  as  1-irregular  meshes  (see  [14]).  The  precise  mathematical  develop¬ 
ment  of  interpolation  spaces  for  1-irregular  meshes  of  bilinear  elements  is  given  in  [14] 
(see  [15]  for  the  biquadratic  elements).  A  systematic  approach  for  the  construction  of 
hierarchic  variable-p  interpolation  spaces  for  1-irregular  meshes  is  given  in  [16].  In  this 
paper  we  will  restrict  ourselves  to  bilinear  (p  =  1)  and  biquadratic  (p  =  2)  elements  of 
the  Lagrangian  family.  In  meshes  of  Lagrangian  elements  vdth  l-irregular  connections 
there  are  two  types  of  nodes,  the  proper  (or  active)  nodes  and  the  improper  (or  con¬ 
strained)  nodes.  For  bilinear  elements  a  proper  node  is  a  vertex  for  all  the  elements 


which  share  the  node  as  shown  in  Fig.  la  (see  also  [14]).  For  biquadratic  elements  a 
proper  node  is  either  a  vertex  for  all  the  elements  which  share  the  node  or  a  midside 
node  for  an  element  and  a  comer  node  for  other  two  elements  which  are  connected  to 
the  first  element  through  that  side  as  shown  in  Fig.  lb  (see  [15]).  The  details  about 
the  implementation  of  the  constrained  approximation  are  discussed  in  [17].  The  con¬ 
strained  nodes  can  be  eliminated  by  introducing  transition  (or  composite)  dements  as 
shown  in  Fig.  Ic  and  Fig.  Id  for  the  meshes  shown  in  Fig.  la  and  Fig.  lb,  respectively. 

We  let  5J,ri,  ==  '^(^)  The  finite  element  solution  it,,  of  the  model 

problem  satisfies:  Find  it,,  6  ^,rjj 

(2.8)  flo(i*fc,Vfc)  =  InM  ^ 

We  let  e,,  :=  u  —  u,,  denote  the  error  in  the  finite  element  solution. 

2.2  Definition  of  f7%-superconvergence  quantities 

Let  T,,  €  T  be  a  finite  element  grid  and  t  £  he  any  element.  Let  F(it)  be 

the  solution-quantity  of  interest,  for  example  we  may  have  F(it)  =  — —  or  — —  or  a 

C/2)  0X2 

linear  combination  of  the  derivatives.  We  now  define  several  geometrical  quantities 
associated  with  the  error  in  the  finite  element  solution  in  the  element  r. 

Given  17,  0  <  q  <  100  we  define  the  following: 

1.  i}%-contour  of  F{y,)  in  the  element  t  for  the  exact  solution  it: 

(2.9)  Cp.^„)(it;r,rfc)  :=  I®  €  r  I  ©(®;F;it,it,„h,r)  =  ^  |, 

Here  @(®;  F;u,u,„/i,r)  is  the  relative  error  as  defined  in  (1.2). 

2.  ri%-hand  of  F(u)  in  the  element  r  €  T,,  for  the  exact  solution  it: 

(2.10)  Bf^^^{u',r,Th)  :=  I*  6  r  I  ©{x;F;u,u^,h,T)  <  } 

3.  Superconvergence  points  of  F{u)  in  the  element  t  E  Tf^  for  the  class  of  exact 
solutions  14: 

(2.11)  r,  Tk)  :=  n  r,  n) 

ueu 
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4-  fi%’Superconveryence  regions  of  F{u)  tn  ike  element  t  €  for  the  class  of  exact 
sobiUons  U: 

(212)  :=  n 

Here  U  denotes  the  set  of  exact  solutions  of  interest. 

K^mark  2.1.  In  some  cases  the  sets  defined  above  may  be  empty. 

Remark  2.2.  The  function  u  is  a  solution  (resp.  U  is  a  class  of  solutions)  of  (2.1) 
for  a  given  set  of  data  (resp.  for  given  classes  of  data).  When  /  =  0  the  solution  u 
satisfies  Laplace’s  equation,  i.e.  —Ait  =  0  (i.e.  the  solution  is  harmonic).  The  major¬ 
ity  of  steady-state  computations  in  engineering  are  done  to  approximate  “harmonic” 
solutions.  (By  “harmonic”  solutions  we  mean  solutions  of  the  homogeneous  differen¬ 
tial  equation  (or  system  of  differential  equations)  with  non-homogeneous  boundary- 
conditions.)  Thus  it  is  important  to  study  superconvergence  for  this  class  of  solutions. 

One  is  interested  to  know  a-priori  the  q%-superconvergence  regions  and  the  super¬ 
convergence  points,  if  they  exist,  for  classes  of  solutions  of  interest  (In  the  practical 
computations  in  plane  elasticity  and  heat-conduction  the  class  of  solutions  of  interest  is 
the  class  of  “harmonic”  solutions  with  a  finite  number  of  algebraic  point  singularities  of 
the  t3rpe  r“).  In  general,  if  7,^  is  any  grid,  it  is  impossible  to  predict  the  locations  of  the 
superconvergence  regions  and  points.  Here  we  will  make  additional  assumptions  about 
the  approximation  which  will  enable  us  to  determine  a-priori  the  q%-superconvergence 
regions  and  points,  if  they  exist.  In  particular  we  will  assume  that: 

1.  We  are  interested  in  the  asymptotic  locations  of  the  i7%-8uperconvergence  regions 
as  a  mesh-parameter  h  tends  to  zero  (details  about  how  the  limit  is  taken  are 
given  below); 

2.  The  global  modes  of  the  error  (pollutions)  are  negligible  when  compared  with 
the  magnitude  of  the  error  in  the  local  best-approximation  (see  Section  3); 

3.  We  are  interested  in  a  specific  class  of  solutions.  The  majority  of  the  results  in 
this  paper  will  be  given  for  the  class  of  harmonic  solutions. 

In  this  paper  we  will  determine  the  asymptotic  q%-superconvergence  regions  and 
points  for  elements  in  the  interior  of  the  mesh  and  smooth  solutions.  (The  solutions 
can  have  point-singularities  at  the  boundaries  but  they  are  analytic  in  the  interior  of 
homogeneous  domains;  all  the  practical  solutions  in  orthotropic  heat-conduction  and 
elasticity  are  of  this  type.)  The  theoretical  residts  are  given  in  [1]  for  a  special  class 
of  locally  periodic  grids.  Here  we  will  demonstrate  through  numerical  examples  that 
the  conclusions  of  the  theoretical  study  hold,  for  all  practical  purposes,  for  the  types 
of  grids  of  quadrilaterals  which  are  used  in  practical  engineering  computations. 
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2.3  The  class  of  locally  periodic  meshes 

We  now  present  the  definition  of  a  special  class  of  locally  periodic  meshes.  Let  ns 
consider  a  locaUy  periodic  grid  defined  as  follows.  Let  0  <  IT  <  JST®  ,  =  (*J,  e  fl , 

(2.13)  =  <=1,2} 

and  assume  that  is  sufficiently  small  such  that  5(0",  H^)  C  fl.  Further,  let  7  be  a 
set  of  multi-indices  (t,  j) ,  =  (*1*“*^  *j €  ft  and 

(2.14)  c(«W),)l):=S(»<«>,4)C5(.”,B), 

be  the  set  of  the  h-cells  (or  cells)  which  cover  exactly  5(0°,  H)  i.e. 

(2.15a)  U  c(«<‘’«,A)  =  5(«»,B) 

(*j)€7 

(2.155)  a(«(‘>,<'),5)nc(«(‘>'*>,5)  =  0  for  (<.,J.)  #  (i„j,) 

We  will  refer  to  5(0°,  IT)  as  the  subdomain  of  periodicity  of  the  mesh  centered  at  0°. 
Denoting  by 

(2.16)  c:=  5(0,1):=  |(*i,ia)|  |ii|  <  1 ,  !®j|  <  l| 

the  unit-  {master-)  cell  c,  every  h-cell  is  an  h-scaled  and  translated  master-cell. 

Let  T  be  a  mesh  on  the  master-cell  (the  master-mesh)  and  be  the  mesh 

on  c{x^*'^\h)  which  is  the  scaled  and  translated  image  of  f.  We  will  consider  the 
ftunily  T  of  locally  periodic  meshes.  Let  Tf^  E  T  and  Tf^{x^,H)  be  the  restriction 
of  on  S{x,H)  and  the  restriction  of  Tf^{x°,H)  on  c(0(’’^),h).  We  assume 
that  ^  (t^  j)  £  -y  i,e.  Tf^{x°,H)  is  made  by  the  periodic  repetition  of  the 

h- scaled  master  mesh. 

The  type  of  meshes  under  consideration  is  depicted  in  Fig.  2a,  where  the  subdo¬ 
mains  5(0°,  H)  D  5(0°,  fTi)  are  covered  by  a  periodic  array  of  cells  shown  in  Fig.  2b 
(where  the  master-mesh  in  the  master-cell  is  shown  in  Fig.  2c).  Outside  the  square 
5(0°,  5^)  the  mesh  is  arbitrary;  it  could  have  curved  elements,  etc.  The  values  of 
depend  on  h  i.e.  f  ^  =  H^{h).  More  precisely  we  will  assume  that  there  exist  constants 
and  C2  independent  of  h  such  that 

(2.17)  Cjfrf  <h<  ,  a  >  1 
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3  The  theoretical  setting 

Given  a  function  tt  and  the  multi-index  a  :=  (ai,  aa)  we  define 


(3.1a) 


:= 


\a\  :=  ttj  -I-  Oj 


(3.16)  (T^u){m)  :=  |D“u|’^(«)j^  ,  As  >  0,  integer 


and  let  us  denote 


(3.2a) 

(3.26) 

1  “  1 ll®lli«(s(«»,jsr^ )) 

(3.2c) 

II  Il5(»«,frj )  •=  53  1  1 

|a|=l 


Let  Q  be  a  polynomial  of  degree  p-f- 1  on  the  master-cell  c  and  let  T  be  the  master- 
mesh.  Then  denote 


(3.3)  P-Q- 

where  is  the  interpolant  of  degree  p  of  the  function  Q  defined  over  the  master- 
mesh  T  (for  which  h  =  1).  As  we  mentioned  above,  any  polynomial  of  degree  p  on  an 
element  Tj,  belongs  to  ^(t*).  Hence  any  polynomial  of  degree  p  on  S{x°,H)  belongs 
to  Si(S(x°,H)).  Hence  p  defined  in  (3.3)  is  c-periodic;  this  can  be  shown  exactly  as 
in  [18].  We  have 


(3.4a) 


^(^>*2)  ~  ^3®a)»  l®jl  1 


(3.46) 


p(xi,l)=  p(xi,-l),  |iil<l 


Further  let 


(3.5)  ^PE»(c)  :=  |tt  6  H^{c)  |  it  satisfies  (3.4)  | 
and 

(3.6)  ^4>br(c)  :=  |«  €  ^per(c)  j  i*U  ^  ‘ST(^)  V  f  6  T  I 
Farther  let 


(3.7o)  =  Bi{p,v)  V  V  €  ^^br(c) 

and 


(3.74)  jf  (,  -  z')  =  0 

Note  that  the  function  exists  and  is  uniquely  determined  (we  will  compute  it  nu¬ 
merically  in  the  examples).  Let  us  also  define  ^  €  H^{c)  by 


(3.8)  V'=  p  —  z^  =  Q  —  w  where  w  :=  +  z^ 

We  will  now  outline  the  main  theorem  of  the  paper.  We  will  make  the  following 
assumptions  about  the  exact  solution  u: 

Assumption  I 

On  S{x°,H) 

(3.9)  I  D°u  I  <  iiT  <  oo,  0  <  |a|  <  p  +  2 

Remark  S.l.  Assumption  I  states  that  the  solution  is  locally  smooth  in  the  subdomain 
S(x°,H)  i.e.  the  subdomain  should  be  sufficiently  far  from  boundaries,  material- 
interfaces  and  points  where  the  data  are  rough. 

Assumption  II 

If  o„  :=  (I>“«)  (®“),  a  =  (ai,a2) ,  |a|  <  p  -f- 1  then 
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(3.10) 


ii>=  ^  al>0 

Farther,  we  assume  that  the  mesh  T(n,h)  is  such  that: 

Assumption  III 

On5(»®,iri),  E^<H<H^ 

(3-11)  l|eklU(»«A)  ^  Ch^Bx 

with  /S  >  (p  +  1)  —  c,  where  c  is  specified  in  the  theorem  below  and  where  C  is 
independent  of  T(f2,  h),  but  it  depends  on  fiT,  H  and  R.  We  will  assume  everywhere 
below  that  H  is  sufficiently  small  depending  on  K,  R  and  p. 

Remark  S.S.  We  do  not  assume  that  u  is  smooth  in  (1  outside  of  5(8$”,  IT).  For 
example,  Q  can  have  a  boimdary  with  reentrant  comers  (as  in  Fig.  2(a))  and  hence  u 
can  be  unsmooth  in  the  neighborhood  of  these  comers.  Nevertheless  assumption  III 
makes  an  implicit  requirement  on  the  (refinement  of  the)  mesh  in  the  neighborhood  of 
these  comers.  If  u  is  smooth  in  a  convex  Q  and  the  mesh  is  quasi-uniform  then 


(3-12)  r>0 

and  hence  in  assumption  III  we  can  take  ^  <  p  +  1  arbitrary. 

Remark  S.S.  Assumptions  II,  III  imply  that  the  principal  part  of  the  error  in  5(8S°,  Hi) 
is  related  to  the  non-zero  (p  +  l)-derivatives  of  the  exact  solution  at  a;'’. 

Let  h))  be  the  function  defined  above,  scaled  and  translated 

^  o  the  cell  of  the  mesh  in  5(a!®,  JJ)  i.e. 

(3.13)  V’k(*):=  ^’^^V'(*),  ^  *  = 

where  x  €  c{x^*’^\h)  and  i  =  i  (as  —  as^***^) . 

h 

It  is  easy  to  see  that  can  be  periodically  extended  over  5(aB“,  Hi).  In  [1]  we  have 
proven  the  following  theorem  based  on  the  theory  of  interior  estimates  (see  [19]-[24]): 

Theorem  1.  Let  Hi  <  H  <  H°  and  the  assumptions  I-III  and  (2.17)  hold  with 

(3.14)  “  =  ^.  ‘'=6(6;rTT)'  '  =  ■' 
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Then  for  any  m  €  S(m^,H^) 

(’•»')  l^w|=l^w|+i*i“'*' 

with  |A|  <  1  and  C  independent  of  h. 

Remark  S.4>  Note  that  ^  (and  thus  depends  on  the  polynomial  Q  which  approx¬ 
imates  well  the  solution  in  5(a‘*,  H).  Theorem  1  states  that  for  small  h  we  have  on 


(3.16)  ^(.)  «  ^(.) 
or  more  precisely 

(3.17)  1^  (.)|  =  I*'  ^  (*)[  +  |A|  Ch^‘  . 

On  the  other  hand  it  is  easy  to  see  that 

(3.18)  1  .) 

'  '  \  •  dx-i  *  «(■<••»•.*)  •  I  e(«<**3)A)  / 

where  C  depends  on  the  constant  R  defined  in  (3.10).  Nevertheless  this  does  not  mean 
that 


(3.19) 

The  case  that 


<=1,2 


(3.20) 


O 

j  I  <  <rQ  >  0 ,  for  dther  t  =  1  or  2 , 


is  very  exceptional.  Hence  we  will  assume  that  (3.20)  does  not  happen  i.e.  we  will 
assume  that  there  exists  a  constant  c  >  0  independent  of  h  such  that 


(3.21) 


dXi 
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This  assumption  could  formally  be  improved  by  imposing  additional  assumptions  on 
the  a^’s  (defined  in  assumption  II).  Note  also  that  (3.21)  is  always  assumed  to  bold  in 
the  classical  superconvergence  theories. 

Remark  S.5.  Under  assumption  (3.21)  Theorem  1  states  that:  A  point  i  tn  the  element 

du 

f  t5  asymptotically  ri%-superconvergence  point  for  F{u)  =  -= — ,  «  =  1,  2,  if  and  only  if 

dXi 

(3.22o)  e(i;/’;Q,t&,l,f)  < 

where 

(3.22h)  Q,w, l,f)  := 

and  zero  otherwise  (see  also  (1.2)). 

Remark  3.6.  Theorem  1  states  that  at  a  superconvergence  point  the  accuracy  is  by  a 
factor  hy  (with  v  given  in  (3.14))  better  than  in  the  other  points.  The  value  of  the 
factor  V  in  (3.14)  is  a  theoretical  one  and  practically  we  usually  see  t/  =  1,  or  higher, 
but  the  theoretical  determination  of  the  largest  v  is  out  of  scope  of  this  paper.  For  some 
examples  of  observed  values  of  v  in  numerical  experiments  see  Section  5.2.  The  case 
(/  >  0  is  of  course  sufficient  for  the  analysis  of  the  existence  of  the  superconvergence 
points  and  17%-superconvergence  regions. 

4  The  methodology  for  determining  the  7/%- 
superconvergence  regions 

4.1  The  method  of  freezing  the  periodicity 

In  this  paper  superconvergence  is  treated  as  local  behavior  and  is  based  on  the 
local  behavior  of  the  solution  in  the  interior  of  the  domain.  We  will  consider  the  class 
of  solutions  which  are  locally  smooth  in  S{x^,H),  namely, 


if  W):=  0 


(4.1)  :=  {11  e  F>(n)|  1 1>“«  <  Jf ,  0<H<p  +  2} 

where  S{x°,H)  denotes  an  interior  subdomain  of  interest  in  which  the  mesh  is  locally 
periodic  as  described  above  (the  subdomain  must  be  a  finite  distance  away  horn  the 
boundary  and  points  of  roughness  of  the  source  term;  see  Fig.  2a).  In  many  instances 
we  are  only  interested  in  the  subclass  of  solutions  in  which  are  harmonic,  namely. 
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(4.2) 


«®  :=  |ii  €  I  At.  =  0  in  ol 

We  may  also  assume  that  the  {unctions  are  harmonic  in  a  subdomain  which  is  slightly 
bigger  than  5(e'’,  H)  and  which  includes  5(c'*,  H)  in  its  interior. 

a.  Tht  classes  of  (p  +  l)-deyree  monomial  aohiUona. 

Let  us  assume  that  for  a  given  locally  periodic  grid  with  corresponding  periodic 
master-mesh  T,  given  material  orthotropy  and  given  class  of  smooth  solutions  li  we 
consider 


(4.3)  Q  :=  {<3  I  S 

fc=l  i=o  ^ 

the  class  of  (p  -I-  l)-degree  monomials  which  occur  in  all  (p  -|-  l)-degree  Taylor-series 
expansions  of  functions  from  U.  Here  k  =  l,...,nd  denotes  a  set  of  linearly 
independent  monomials  which  form  a  basis  for  Q.  For  example,  let  us  assume  that  U 
is  the  class  of  smooth  solutions  given  in  (4.1);  in  this  case  we  may  choose 

(4.4)  (?;^(xi,®2)  :=  ,  l<h<nd  =  p-|-2 

and  we  obtain  the  class  of  all  (p  -h  l)-degree  monomials  Q^. 

In  the  case  that  we  are  interested  only  in  the  class  of  harmonic  solutions  14^  given 
in  (4.2)  we  will  take  Q  as  the  two-dimensional  linear  space  of  harmonic  monomials  of 
degree  (p  -|- 1)  denoted  by  namely, 


(4.5o)  I  “fc  }  > 

h=l  * 

(4.56)  Qf  (®i,a!a)  =  ,  Qf(®i,x,)  =  Im(z*^^) ,  z  =  x^  -|-  tx,  . 


In  the  previous  Section  we  outlined  Theorem  1  which  states  that  we  can  obtain 
the  asymptotic  values  of  the  error  for  any  smooth  solution  u  in  the  interior  of  a  peri¬ 
odic  mesh- subdomain  by  solving  the  periodic  boundary-value  problem  (3.7),  using  the 
master-mesh  T  over  the  master-cell  c,  with  data  obtained  from  the  local  (p-|-  l)-degree 
Taylor-series  expansion  Q  of  the  exact  solution.  Based  on  this  result  we  will  construct  a 
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numerical  procedure  to  determine  the  i7%-8uperconTergence  quantities  for  a  given  class 
of  smooth  solutions  by  employing  the  corresponding  class  of  (p+  l)-degree  monomials. 

In  order  to  apply  the  results  of  the  theoretical  study  to  the  practical  meshes,  for 
which  the  mesh  is  not  locally  periodic  (like  for  example  the  mesh  shown  in  Fig.  3a), 
the  following  technique  of  freezing  the  periodicity  will  be  employed: 

1.  Let  Uq  be  an  interior  patch  of  elements  of  interest  (shown  shaded  gray  in  Fig.  3a) 
in  which  we  study  the  q%-superconvergence.  Define  the  s-layered  patch  of  ele¬ 
ments  s  >  1,  surrounding  the  patch  Uq  by 

(4.6)  :=  U  “'x  1  :=  U 

X€JV(t')  xeJV(T') 

where  N{t*)  denotes  the  set  of  the  vertices  of  element  r',  0;^^  is  the  patch  of 
the  elements  connected  to  vertex  X.  The  patch  is  shown  with  thick  black 
perigram  in  Fig.  3b. 

2.  Complete  the  patch  to  a  periodic-grid  over  a  slightly  larger  square  periodic¬ 
cell  which  encloses  the  patch  as  shown  in  Fig.  3c.  The  periodic-cell  is  then  scaled 
and  translated  to  the  unit  master-cell  c. 

3.  Assume  that  the  mesh  in  the  neighborhood  of  element  r  is  made  from  the  periodic 
repetition  of  h-cells  obtained  from  the  master-cell  (by  h-scaling  and  translation) 
and  let  h  tend  to  zero.  Then  the  theoretical-setting  of  the  previous  Section 
applies  and  the  asymptotic  error-function  in  the  mesh- cell  (Oq  can  be  obtained  by 
solving  the  periodic  boundary- value  problem  over  the  master-cell  c.  Based  on  the 
results  of  the  theoretical  study,  outlined  in  Section  3,  the  asymptotic  locations 
of  the  17%-superconvergence  quantities  in  any  element  t  C  for  a  solution  u 
(resp.  class  of  solutions  U)  can  be  obtained  from  the  corresponding  quantities 
(defined  on  the  master-cell  c)  for  the  local  (p-H  l)-degree  Taylor-series  expansion 
Q  of  u  (resp.  the  class  Q  of  the  local  (p  -f-  l)-degtee  Taylor-series  expansions  of 
functions  u  £  14),  namely: 

(4.7)  Um 


(4-8)  Um  ej* ,(»;  r,  TJ  =  BJ*  ,(Q:  f ,  f ) 
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(4.9) 

(4.10) 

Bm  *J5,,(W;T,rt)=  Rj*)(C;f,T) 

Remark  The  above  limits  hold  for  the  locally  periodic  meshes  under  the  assump¬ 
tions  of  the  theoretical  analysis  outlined  in  Section  3  (see  [1]  for  the  details).  Hence 
for  the  general  grids  the  limit  should  be  understood  for  the  mesh  which  is  constructed 
by  freezing  the  periodicity  (as  is,  for  example,  shown  in  Fig.  3c). 

In  the  numerical  examples  it  will  be  demonstrated  that  the  results  of  i7%-superconver- 
gence  obtained  from  the  above  methodology  hold,  for  all  practical  purposes,  for  the 
complex  grids  used  in  en^eering  computations  (provided  that  the  pollution-error  is 
controlled,  the  approximation  is  in  the  asymptotic  range  and  that  sufficient  number  of 
mesh-layers  are  included  in  the  patch  w^). 

Remark  4.2.  Above  we  assumed  that  Uq  =  of  course  we  can  use  smy  other  patch 
as  {Jq,  for  example  consisting  of  one  element  only. 

4.2  Determination  of  the  i7%*superconvergence  quantities 

The  asjrmptotic  i7%-contours  for  a  given  solution  u  can  be  obtained  by  contouring  the 
function  defined  in  (3.8),  corresponding  to  the  local  Taylor-series  expansion  Q  of 
the  solution  u.  The  superconvergence  points  i  for  a  given  class  of  solutions  U  satisfy 


(4.11)  F(f/>iXi)=  0,  l<t<nd 

i.e.  when  there  is  a  superconvergence  point  then  the  zero-contours  of  F(^i )  intersect 
at  i  for  1  <  t  <  nd.  Here  V’i  '=  Pi  —  which  is  obtained  from  (3.7)  for  = 

—  (Qj)?*^  where  is  the  i-th  basis  monomial  of  the  nd-dimensional  monomial 

space  Q  corresponding  to  the  class  as  discussed  in  Section  4.1.  We  also  let  := 

The  as3rmptotic  f7%-superconvergence  bands  for  a  solution  u  can  be  determined 
from  the  function  F'(V’)  by  using  piecewise  linear  interpolation  of  F{^)  on  a  sufficiently 
refined  uniform  mesh  obtained  by  subdividing  the  quadrilateral  of  interest. 

The  asymptotic  q%-superconvergence  regions  for  a  class  of  solutions  U  can  be  de¬ 
termined  by  using  numerical  optimization.  In  particular,  consider  the  uniform  subdivi¬ 
sion  of  the  element  t  into  subtriangles  with  vertices  at  the  set  of  points  E  :=  • 

Define  the  relative-error  function  at  the  point 
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(4.12)  F\  Q,  {wjfei,  1,^)  :=  max - *=^5 - ’ 

•  m«  |  Ea,n^i)(^,)| 

j— ii—i")*  t=i 

Then  the  function  6i{'^i}{^itlf‘r)  can  be  defined  for  any  point  e  E  f  by 

using  linear  interpolation  in  the  subtriangles.  The  asymptotic  i7%-superconvergence 
regions  in  the  element  r  can  be  approximated  using  the  level-sets  of  the  functions 
®l(u)(«;  F\  Q,  l,f )  i.c. 

(4.13)  )(fi;‘r,T)  «  |«  €  f  I  ©!(„)(», •^;  G,  l,f)  <  i7%| 

We  will  call  the  above  approach  the  dirtci  approach.  It  is  also  possible  to  use  a  simplified 
approach  which  avoids  the  use  of  numerical  optimization  at  every  point.  First  observe 
that 


max 


I  S  aiF(V’iX*)l  ' 


»=l 


(4.14) 


m»x  |  ra(m)«j)l 


min 


1=1 


/J 


-1 


Hence  we  can  define  the  function 


(4.15) 
where 

(4.16) 


)(«))’. 


1=1 


: 


=  imn 

a.. 


'"  .max  I  EaiWi)(^j)l'^ 


V 


nd  , 

(E°?)> 

1=1 
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The  quantity  Zs  can  be  computed  using  numerical  optimisation.  Let 

(4.17)  :=  {i  €  r  I  } 

denote  the  approximate  regions  of  q%*superconvergence  for  the  class  of  solutions  Q 
obtained  by  the  simplified  approach. 

The  common  q%-superconTergence  region  in  an  element  r  for  a  given  class  of  meshes 
and  materials  can  be  determined  by  using  either  the  direct  or  the  simplified  approach. 

Remark  4-3.  Note  that  we  have 

(4.18) 

Therefore  the  simplified  approach  results  to  a  conservative  estimate  for  the  i7%-superconvergence 
regions.  For  the  class  of  harmonic  functions  we  observed  that  the  regions  obtained 
by  the  two  approaches  are  very  similar. 

Remark  4-4-  The  functions  defined  in  (4.12),  (4.15)  depend  on  the  set  of  points  H. 

To  ensure  good  accuracy  in  the  approximation  of  the  17%-superconvergence  regions  a 
sufficient  number  of  points  must  be  employed.  In  the  examples  below  the  set  of  points 
H  was  obtained  by  subdividing  the  quadrilateral  into  two  triangles  and  by  refining  these 
triangles  uniformly  several  times.  The  points  in  H  are  the  vertices  of  this  triangulation. 

5  Numerical  studies  of  super  convergence 

We  win  now  give  examples  of  numerical  studies  of  classical  superconvergence 
and  1/%-superconvergence  for  complex  grids  of  quadrilaterals.  In  this  Section  we  will 
address  the  following  questions: 

1.  Given  any  periodic  mesh  of  quadrilaterals  (with  complex  geometry  and/or  local 
refinements  etc.),  where  are  the  superconvergence  points  for  the  derivatives,  for 
the  class  of  harmonic  and  general  solutions? 

2.  Given  an  element  in  the  interior  of  a  complex  grid  of  quadrilaterals,  where  are 
the  points  of  least-error  for  the  derivatives  for  the  class  of  harmonic  and  general 
solutions? 

3.  Can  we  put  some  reasonable  restrictions  on  the  grid  (which  can  be  controlled  in 
the  computation)  such  that  q%-8uperconvergence  for  the  derivatives  is  assured 
for  any  solution  (including  solutions  with  algebraic  point-singularities  i.e.  all  the 
solutions  which  are  of  interest  in  the  practical  computations)  in  the  interior  of 
any  practical  grid? 
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5.1  Determination  of  the  superconvergence  points  for  the 
periodic  meshes 

In  [1]  we  demonstrated  that  the  classical  snperconvergence  points  for  the  derivatives 
in  periodic  meshes  of  triai^es  or  squares  can  be  determined  numerically.  Here  we 
use  the  same  methodology  to  find  the  superconvergence  points  for  periodic  meshes  of 
square  elements  with  refinements  (or  meshes  obtained  from  these  meshes  via  a  smooth 
transformation).  We  studied  in  detail  the  superconvergence  in  the  periodic  meshes 
shown  in  Fig.  4  and  we  give  the  results  below. 

5.1.1  Grid  with  mesh-interface  with  constrained  nodes 

In  order  to  study  the  superconvergence  in  the  elements  near  a  mesh-interface  with 
constrained  nodes  we  considered  the  periodic-cdl  shown  in  Fig.  4a.  Note  that  although 
the  theoretical  setting  was  outlined  for  a  square  periodic  unit-cell  c,  it  is  also  valid  for 
the  rectangular  unit-cell  employed  in  this  example.  Here  c  =  (— 0.5h,0.5h)x  (0,1)  and 
h  =  0.05.  We  are  interested  in  determining  the  superconvergence  points  for  the  x^-  and 
Xj-derivatives  for  the  class  of  harmonic  and  the  class  of  general  solutions  for  bilinear 
and  biquadratic  elements. 

(a)  Bilinear  elements 

For  bilinear  elements,  the  only  quadratic  harmonic  monomial  which  gives  non-zero 
error  is  <?f(®i,*2)  =  *?  —  ®j*  Therefore  we  have  superconvergence  of  the  Xj-  and  x,- 

a.,  ’  a*, 

error-function  which  corresponds  to  )  which  are  shown  in  Fig.  5a.  The  points  of 
intersection  of  these  lines  are  superconvergence  points  for  the  x^-  and  Xj-derivatives, 
simultaneously;  these  points  are  given  in  Table  1. 

Further  we  note  that  for  the  class  of  general  solutions  there  are  two  general  quadratic 
monomials  which  give  non-zero  error,  namely  Qf  (sBi,  SB2)  ~  >  ^2)  =  ^2* 

denote  the  error-functions  in  the  periodic  mesh-cell,  corresponding  to  these  monomials, 
by  and  respectively  and  analyze  their  structure.  First  consider  the  function 
and  let  us  denote  by  p,  the  interpolation-error  function  introduced  in  (3.3)  which 
is  associated  to  Qp.  We  can  easily  see  that  Qp  and  (Qp)P^^  are  independent  of  x^ 
and  that  p,  =  0  on  every  boundary  line  of  the  elements  parallel  to  the  mesh-interface. 
In  addition  for  every  fixed  Xj,  (<?p)f*^  is  a  linear  interpolant  of  Qp.  From  this  we  see 
that  in  (3.7a)  we  have 

(5.1)  ^c(P2.v)  =  0  V  v  €  5J,per(c) 

and  hence  is  a  constant  function.  Therefore 


derivatives  along  the  zero-lines 


,  respectively  (where  V'P  is  the  periodic 
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(5.2) 


»  =  1,2 


dVf  _  ^ 

dZi  dxi  * 


and  (as  special  case)  vanishes  identically. 

asj 

Let  ns  now  analyse  Because  is  symmetric  with  inspect  to  Xj,  so  is  s'>. 
Hence  in  every  element  of  sixe  h,  is  constant  for  fixed  x,.  In  addition  it  is  easy 
to  write  difference  equations  for  the  values  of  at  the  nodal  points  (of  dements  of 
size  h).  From  this  we  condude  that  the  function  has  to  be  linear  in  the  subdomain 
occupied  by  dements  of  size  h  and  using  the  periodidty  we  condude  that  x^i  has  to 
be  constant  in  this  subdomain.  From  the  difference  equations  it  is  easy  to  see  that  the 

function  is  not  constant  in  the  subdomain  occupied  by  dements  of  size  — .  A  more 

detailed  analysis  shows  that  the  derivatives  of  the  function  decay  exponentially 
with  the  distance  from  the  mesh-in^erface. 

Because  is  constant  and  vanishes  identically  as  discussed  above,  there  are 

OX\ 

superconvergence  lines  for  the  Xj*derivatives;  these  lines  are  the  same  as  the  supercon¬ 
vergence  lines  for  the  x^-derivatives  for  the  harmonic  solutions  (shown  by  thick  line  in 
Fig.  5a).  For  the  x^-derivative  there  is  one  superconvergence  point  in  each  dement  at 

the  intersections  of  the  zero-lines  of  -5-^,  -5-^,  as  shown  in  Fig.  5b. 

Since  2^< ,  t  =  1,  2,  is  constant  in  the  dements  of  size  h  the  simultaneous  super¬ 
convergence  points  for  both  derivatives  are  located  at  the  centers  of  these  dements. 

We  note  that  in  the  dements  of  size  —  the  simultaneous  superconvergence  points  are 
not  at  the  center  of  the  dement  (for  harmonic  solutions)  or  they  do  not  exist  at  all 
(for  general  solutions).  Neverthdess  in  the  dement  of  size  —  not  at  the  interlace  the 

q%- value  at  the  centers  of  the  dements  is  very  small  and  these  points  are  simultaneous 
superconvergence  points  for  very  small  tolerance.  This  follows  from  the  exponential 
decay  of  the  derivatives  of  the  function  z'>  with  the  distance  from  the  interface. 


(b)  Biquadratic  elements 

For  the  biquadratic  dements  we  note  that  the  harmonic  and  general  monomials  of 
degree  3  (namely  (x^, Xj)  =  *?  3xjx|,  Qp(3!i, s^s)  ~  *5  ~  and  Qp(z^,  Xj)  = 

—  ^\)  result  in  identical  error-functions  V'P)  ill  dl  the  dements. 

Let  us  analyze  the  functions  V'Pi  *  =  li  2,  analogously  as  in  the  case  of  bilinear 
dements.  First,  let  us  address  the  case  t  =  2.  As  in  the  case  of  bilinear  dements 
Q2  and  (Qp)f^^  are  independent  of  x^  and  p,  is  identically  zero  on  the  dement- 
boundaries  and  the  center-lines  paralld  to  the  mesh-interface.  Moreover  (Qp)P^^  is  a 
quadratic  interpolant  of  Qp  for  every  fixed  x^.  Taking  into  account  that  = 

X],  analogously  as  before,  we  condude  that 
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(5.3)  B^{p^fV)  —  0  V  t>  € 

and  hence  is  constant.  Therefore 


.*-12 

dxi~  dxi^  ’ 


and  vanishes  identically  and  hence  we  have  snperconvergence  lines  for  the  x.- 

OXi 

derivative  which  are  the  zero-lines  of  (these  lines  are  the  thick  lines  shown  in 

oxi 

Fig.  5c).  Farther  the  zero-lines  of  -r —  are  the  Gauss-lines  parallel  to  the  mesh- 

oX] 

interface. 

Let  us  now  consider  the  case  t  =  1.  First  we  observe  that  is  antisymmetric  with 
respect  to  x^  and  hence  z^>(0,x,)  =  0.  In  the  elements  of  size  h  the  antisymmetry  of 
z^\  implies  that  z'l  has  to  be  linear  in  x^.  From  the  periodicity  of  it  follows  that 

is  identically  zero.  Therefore,  in  the  elements  of  size  h,  is  identically  zero  and 

OX] 

the  superconvergence-lines  for  the  x^-derivative  in  the  elements  of  size  h  are  the  Gauss- 
lines  parallel  to  the  mesh-interface  (the  zero-lines  of  -s-^).  The  superconvergence  lines 

C/X] 

for  the  x^-derivative  in  the  elements  of  size  h  are  the  Gauss-lines  perpendicular  to  the 
mesh-interface  (the  zero-lines  of  -r-^). 

OXi 

For  the  elements  of  size  —  the  situation  is  more  complicated  because  is  non- 

zero  although  its  derivatives  decay  exponentially  with  the  distance  from  the  mesh- 
interface.  Hence  the  superconvergence  lines  for  the  Xj -derivative  in  the  elements  of 
h  d^? 

size  —  (the  zero-lines  of  -5 — )  are  curved  as  shown  in  Fig.  5c.  For  elements  which  are 

2  ijX\ 

far  from  the  mesh-interface  the  superconvergence-lines  are  practically  the  Gauss-lines 

9il>?  h 

perpendicular  to  the  mesh-interface.  Further  — —  ^  0,  in  the  elements  of  size  — ,  and 

UX3  2 

the  values  of  17%  at  the  Gauss-lines  parallel  to  the  mesh-interface  (the  zero-lines  of 

—  )  do  not  exceed  0.001%.  Therefore  the  Gauss-lines  parallel  to  the  mesh-interface 
dx2 

are  17%-superconvergence  lines  for  the  x, -derivative  with  tolerance  less  than  0.001%. 

From  what  we  said  above  we  see  that  there  exist  “almost”  simultaneous  supercon¬ 
vergence  points  (with  tolerance  less  than  0.001%)  and  hence  for  all  practical  purposes 
we  can  use  them  as  superconvergence  points  (these  points  are  listed  in  Table  1).  This 
was  also  mentioned  in  Remark  1.1. 
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5.1.2  Grid  with  mesh-interface  with  transition  elements 


We  also  studied  the  superconvergence  in  the  elements  near  a  mesh-interface  with 
transition  elements  by  employing  the  rectangular  periodic  mesh-cell  shown  in  Fig.  4b. 

(a)  Linear  elements 

In  order  to  find  the  superconvergence  points  for  the  class  of  harmonic  solutions 
we  considered  the  periodic  error  functions  ifif  which  correspond  to  the  quadratic 
harmonic  basis  monomials,  =  af  —  al,  function 

z^2  is  antis3rmmetric  with  respect  to  the  axis  =  0  (note  that  the  mesh  is  symmetric 
about  this  axis).  Using  the  periodicity  condition  we  see  that  2^3  =  0.  Hence  is 
identically  equal  to  zero  at  the  nodes  and  thus  in  the  bilinear  elements  there  exist 

lines  of  superconvergence  for  the  and  ^^-derivatives  at  the  zero-lines  of  — — ,  ■  , 

OX\  UX3 

respectively.  In  the  triangular  elements  there  exist  superconvergence  points  for  the  x^- 

8il>? 

(resp.  X,-)  derivative  at  the  intersection  of  zero-lines  of  -v--',  -r —  (resp.  ,  — — ); 

C/X^  vX^  vXj 

these  points  are  given  in  Table  2a.  For  the  class  of  general  solutions  there  are  three  basis 
monomials  of  quadratic  degree,  Qf(x^,x,)  =  xf,  —  >31  Qs 

We  determined  the  superconvergence  points  in  the  bilinear  square-elements  for  the  x^- 

diff? 

(resp.  X--)  derivative  at  the  intersection  of  the  zero-lines  of  -3-^,  -3 —  (resp.  -3-^, 
*  axi  axi  0x2 

-3-2-);  these  points  are  given  in  Table  2a.  We  note  that  =  P3  =  0  in  the  square 
0X2 

elements;  in  the  triangular  elements  we  utilized  also  the  function  V'P-  Hi  tlus  Table 
we  also  indicate  which  ones  of  the  superconvergence  points  for  the  class  of  harmonic 
solutions  in  the  triangular  elements  are  also  valid  for  the  class  of  general  solutions. 


(h)  Quadratic  elements 

We  computed  the  superconvergence  points  for  the  Xj-  (resp.  x,-)  derivative  in  the 
elements  for  the  class  of  harmonic  solutions  by  considering  the  intersections  of  the 
zero-lines  of  the  x^-  (resp.  x,-)  derivatives  for  the  error-functions  corresponding  to  the 
two  cubic  harmonic  basis  monomials;  the  points  are  given  in  Table  2b.  Note  that  for 
the  class  of  general  solutions  there  are  no  superconvergence  points  in  any  element  of 
this  mesh. 


5.1.3  A  periodic-mesh  with  seversd  refinements 

We  also  studied  the  superconvergence  for  the  periodic  mesh  of  bilinear  elements 
shown  in  Fig.  4c.  For  this  mesh  we  determined  the  points  of  simultaneous  super- 
convergence  for  the  x^-  and  x^-derivatives  for  the  class  of  harmonic  solutions;  these 
points  are  given  in  Table  3.  Note  that  for  some  elements  with  constrained  nodes  the 
superconvergence  points  for  the  derivatives  coincide  with  a  vertex  of  the  element. 
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5.2  Rate  of  convergence  at  the  superconvergence  points 

We  performed  model  computations  on  finite-element  meshes  with  mesh-interface 
(with  constrained  nodes  or  with  transition  elements)  and  we  checked  the  rate  of  a>n- 
vergence  of  the  derivatives  at  the  superconvergence  points  which  are  given  in  Tables 
1,  2  for  these  meshes.  We  found  that  the  rate  of  convergence  at  the  superconvergence 
points  (given  in  Tables  1,  2)  is  of  order  (p  -f- 1).  We  also  checked  the  rate  of  conver¬ 
gence  at  the  p  X  p  Gauss- Legendre  points  (for  p  =  1,  2)  in  the  square  elements  near 
the  mesh-interface  and  found  that  these  points  are  not  superconvergence  points,  in 
general. 

In  the  computations  we  considered  the  Dirichlet  boundary-value  problem  for  the 
Lapladan  in  the  domain  Q  =  (0, 1)^  with  Dirichlet-data  obtained  from  the  harmonic 
exact  solution  v.{x^,x^)  =  sin(xxj)sinh(wx,).  The  domain  II  was  meshed  with  the  pe¬ 
riodic  grids  with  ceU-size  h  as  shown  in  Figs.  6  and  7.  For  these  meshes  we  computed 

I  dc  I  I  dc  I 

the  maximum  values  of  (the  maximum  was  computed  over  the  supercon¬ 


vergence  points  or  the  Gauss- Legendre  points)  in  the  square-elements  of  size  —  at  the 

mesh-interface  and  one-layer  away  from  it. 

The  results  of  the  convergence  study  for  the  mesh-interface  with  constrained  nodes 
are  given  in  Tables  4a,  4b,  5a,  5b.  In  Table  4a  (resp.  Table  5a)  we  give  the  values  of 

dc  I 

'^  )  at  the  superconvergence  point  and  the  Gauss- Legendre  point  in  the 


^1 


\dxi\ 


elements  of  size  —  at  the  interface  and  one-layer  away  from  the  interface  for  bilineair 
approximation  (p  =  1).  In  Table  4b  (resp.  Table  5b)  we  give  the  maximum  valuer  of 


1^6/, 

>dxi 


(resp.  j^^|)  (computed  over  the  superconvergence  points  or  the  Gauss-Legendre 
points)  in  the  same  elements  for  biquadratic  approximation  (p  =  2).  We  observed  that: 

(i)  When  the  superconvergence  points  are  employed  we  have 


where  denotes  a  superconvergence  point. 

(ii)  When  the  Gauss-Legendre  points  are  employed  we  have 


(5.56) 


I  I  :=  m«  «  Ch- 

'  G—L  I  dxi  ^ 
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(5.5e) 


1^1 

I  a*,  I G-L  • 


max 

mO-t-er 


wheie  denotes  a  Gauss- Legendre  point. 


We  note  that  the  values  of 


dZi 


apprear  to  be  snperconvergent  at  the  Ganss- 


Legendre  points.  (More  precisely  they  are  not  snperconvergent  but  it  cannot  be  seen 
in  the  computations  because  the  points  are  17%-superconvergence  with  very  small  tol¬ 
erance,  as  explained  in  Section  5.1.1.)  We  also  note  that  the  values  of  the  constant 
C  at  the  elements  one-layer  away  from  the  mesh-interface  are  much  smaller  than  the 
values  of  C  in  the  element  at  the  mesh-interface. 

The  results  of  the  convergence  study  for  the  mesh-interface  with  transition  elements 
are  given  in  Tables  6a,  6b,  7.  In  Table  6a  (resp.  6b)  we  give  the  maximum  values  of 

idzi 


(the  maximum  is  computed  either  over  the  snperconvergence  points  or  the  Gauss- 


Legendre  points  in  each  element)  for  linear  (resp.  quadratic)  elements  of  size  —  at  the 
mesh-interface  or  one-layer  away  from  it.  We  observed  that  when  the  snperconvergence 
points  from  Tables  2a,  2b  (resp.  the  Gauss-Legendre  points)  are  employed  the  rate 

I  ^^6 

dz\ 


computed  at  the  superconvergence  points  (given  in  Table  2a)  and  at  the  center  of  the 
element  for  the  triangular  elements  ,  r,  (shown  in  Fig.  4b).  We  observe  that  the  rate 
of  convergence  of  the  z^-derivative  is  (p  -i- 1)  (resp.  p)  at  the  superconvergence  point 
(resp.  at  the  center  of  the  element). 

In  theoretical  studies  of  superconvergence  [4,  5]  the  superconvergence  is  analyzed 
for  discrete  mean-square  norms  evaluated  over  subdomains  (by  employing  the  set  of  all 
snperconvergence  points  in  a  subdomain).  A  simple  generalization  of  the  discrete-norm 
used  in  [4,  5]  and  elsewhere  is 


where  N  denotes  the  number  of  points  {z,}^i  (used  in  the  definition  of  E^)  in  the 
subdomain  for  which  the  discrete-norm  is  evaluated  (for  q  =  2  the  norm  of  [4,  5]  is 
obtained).  For  the  Dirichlet  problem  with  exact  solution  u{x^,X2)  =  sin(irzi)  sinh(7rz2) 

in  the  domain  (2  =  (0, 1)’  and  with  mesh-interface  along  the  line  z,  =  -  (as  shown 

in  Fig.  6a)  we  computed  the  discrete-norm  using  the  Gauss-Legendre  points  in 
the  elements  in  the  subdomain  Gg  =  (0.25,0.75)  x  (0.50,0.75)  for  meshes  of  bilinear 
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ftnd  biquadratic  elements  and  several  values  of  q  £  [1,  oo].  The  values  of  the  discrete- 
norm  and  its  convergence  rates  are  given  in  Tables  8a  (resp.  8b)  for  bilinear  (reap, 
biquadratic)  elements.  We  observe  that 

(5.7)  «  Chr^l 

(See  also  Remark  1.1.) 

Remark  5.1.  As  noted  in  Remark  3.3,  by  emplo3ring  the  theory  (which  is  given  in  detail 
in  [1])  we  cannot  show  that  the  gain  of  the  convergence  rate  is  A"  with  t/  =  1,  as  has 
been  observed  in  the  computations. 

5.3  Relative  error  at  the  superconvergence  points  for  meshes 
which  are  not  locally  periodic. 

The  classical  superconvergence  points  are  obtained  assuming  that  the  grid  is 
locally  periodic.  However,  the  meshes  which  are  used  in  practical  computations  are 
not  expected  to  satisfy  this  assumption.  In  such  meshes  we  can  still  use  the  points 
from  a  periodic  pattern  (for  example  the  points  for  the  mesh-interface  given  above)  to 
sample  the  solution,  if  that  pattern  appears  in  the  interior  of  the  grid.  We  now  show 
that  by  using  such  points  we  obtain  more  accurate  values  of  the  derivatives  than  if  we 
employ  the  p  x  p  Gauss- Legendre  points  in  each  element  (as  is  customarily  done). 

We  considered  the  domain  fl  =  (0, 1)^  with  Neumann  boundary  conditions  g  on  dQ 
and  source-function  / corresponding  to  the  exact  solution,  u(xj,X3)  = 

Here  we  let  A  =  10  smd  (a,  b)  =  (0,0).  A  mesh  of  bilinear  elements  was  generated  using 
an  adaptive  algorithm  as  shown  in  Fig.  8a.  We  will  focus  on  the  elements  in  Patch 
I  (shown  in  Fig.  8a  and  Fig.  8b)  and  Patch  II  (shown  in  Fig.  8c  and  Fig.  8d).  In 

_ 

the  elements  shown  in  each  patch  we  computed  the  relative  error  ©(»;  -5 — ;  rt,  tt^.  A,  r), 

ox^ 

t  =  1,  2,  in  the  x^-  and  Xj-deiivatives  of  at  the  center  of  the  element  (the  Gauss- 
Legendre  point)  or  at  the  superconvergence  point  which  was  obtained  by  assuming 
that  the  element  belongs  to  the  periodic  mesh-interface  pattern  (shown  in  Fig.  4a).  In 
Table  9a  (resp.  Table  9b)  we  give  the  values  of  the  relative  error  in  the  derivatives  at 
the  superconvergence  points  and  the  Gauss-Legendre  points  for  the  elements  in  Patch 
I  (resp.  Patch  II).  We  observe  that  in  the  elements  with  constrained  nodes  in  Patch 
I  (elements  1-16)  the  relative  error  at  the  superconvergence  points  is  substantially 
smaller  than  the  relative  error  at  the  Gauss-Legendre  points.  In  particular  the  relative 
errors  in  the  Xj -derivative  at  the  Gauss-Legendre  points  are  around  25%  whereas  at 
the  superconvergence  points  they  are  around  2%.  In  the  case  of  the  elements  with  con¬ 
strained  nodes  in  Patch  II  (elements  1-4)  the  relative  error  in  the  Xj-derivative  at  the 
superconvergence  points  is  much  smaller  than  the  relative  error  in  the  Xj-derivative  at 
the  Gauss-Legendre  points.  Hence  we  can  conclude  that  the  superconvergence  points 
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from  a  periodic  mesh'pattern  should  be  used  instead  of  the  Gauss* Legendre  points  to 
sample  the  derivatives  of  the  solution  in  the  interior  of  any  general  grid  which  (may 
not  be  locally  periodic  but)  includes  that  pattern  in  its  interior. 

5.4  77%-superconvergeiice  in  periodic  meshes  with  local  re¬ 
finements 

We  shall  now  describe  numerical  examples  to  illustrate  the  methodology  described 
in  Section  4  to  determine  regions  of  17%-superconvergence.  We  shall  first  use  the 
methodology  to  find  i7%-superconvergen(»  regions  t  =  1,  2,  for  the 
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derivatives  in  the  element  r  in  the  periodic  mesh  T  for  the  class  Q  of  (harmonic  or 
general)  polynomial  solutions  of  degree  (p  1).  We  will  then  demonstrate  that  the 
methodology  can  also  be  used  to  find  i}%*superconvergence  regions  in  elements  in  the 
interior  of  any  locally  refined  mesh  or  a  mesh  with  distorted  elements. 

We  considered  the  periodic  meshes  shown  in  Figs.  4a-4c.  In  Figs.  9a-9d  we  give  the 
17%-superconvergence  regions  for  the  elements  in  the  mesh  shown  in  Fig.  4a.  Figs.  9a 
and  9b  show  T),  for  =  5%,  10%,  25%,  for  p  =  1  and  p  =  2,  respectively. 

In  Figs.  9c  and  9d  the  regions  for  i/%  =  5%,  10%,  25%,  are  given  for 
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p  =  1  and  p  =  2,  respectively.  It  can  be  observed  that  the  i7%-superconvergence 
regions  for  the  -derivative  in  the  elements  with  constrained  nodes  are  very  different 
from  the  »7%-superconvergence  regions  for  the  z^-derivative  in  the  elements  away  from 
the  interface.  However,  the  i7%-superconvergence  regions  for  the  Zj-derivative  in  the 
elements  with  constrained  nodes  are  very  similar  to  the  i/%-superconvergence  regions 
for  the  Zj-derivative  in  the  elements  away  from  the  mesh-interface. 

The  1/%-superconvergence  regions  for  the  elements  in  the  mesh  shown  in  Fig.  4b 
are  given  in  Figs.  10a- 10c,  for  the  class  of  harmonic  solutions.  It  can  be  observed 

that  the  effect  of  transition  elements  in  the  elements  of  size  —  decays  after  one  layer 

of  elements.  In  Figs,  lla-lld  we  give  »7%-superconvergence  regions  for  the  class  of 
general  solutions.  Note  that  in  the  triangular  elements,  there  are  no  regions  for  either 
derivative  for  f/%  <  25%  and  that  the  effect  of  the  transition  elements  is  significant  up 
to  4  layers  of  elements  in  the  refined  mesh. 

We  also  studied  the  17%-superconvergence  regions  in  the  elements  in  the  periodic 
mesh-pattem  of  Fig.  4c.  In  Figs.  12a-12d  we  show  the  i}%  superconvergence  regions 
for  p  =  1  and  p  =  2  for  the  class  of  harmonic  solutions.  Figs.  12a  and  12b  show  the 
regions  'RP'^  (G^;t,  7),),  for  ij%  =  5%,  10%,  25%  for  p  =  1  and  p  =  2,  respectively.  In 

Figs.  12c  and  12d  we  show  the  regions  (Q®;t,T^),  for  77%  =  5%,  10%,  25%  for 

p  =  1  and  p  =  2,  respectively.  It  can  be  seen  from  Figs,  12a  and  12c  that  the  values  of 
the  7f%  function  for  the  z^  and  Zj  derivatives  at  the  center  of  the  elements  which  have 
constrained  nodes  is  greater  than  25%. 
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To  demonstrate  the  decay  of  the  effect  of  mesh-interface  on  the  superconvergence 
we  computed  the  values  of  maximum  relative  error  for  the  class  of  harmonic  solutions 
at  the  Gauss- Legendre  points.  In  Table  10a  (resp.  10b)  we  give  the  values  of 


(5.8) 


4®  ^:=  max©(i®  .|!i;Q,,D,l,f) 


Q€fl" 


for  the  elements  f  of  size  — ,  which  are  separated  by  0, 1,  2  and  3  layers  from  the  mesh- 

£t 

interface  for  bilinear  (resp.  biquadratic)  approximation  for  the  mesh-interface  with 
constrained  nodes.  In  Tables  10c  (resp.  lOd)  we  give  the  values  of  for  linear 

(resp.  quadratic)  elements  for  the  mesh-interface  with  transition  elements.  In  both 

h 

cases  we  computed  the  values  only  for  the  square  elements  of  size  — .  We  observed  that 
the  value  of  decays  rapidly  as  we  consider  elements  away  from  the  interface. 


Hence  the  Gauss-Legendre  points  in  the  elements  of  size  —  which  are  separated  by 

several  mesh-layers  from  the  mesh-interface  are  17%-superconvergence  points  for  very 
small  tolerance. 


5.5  77%  superconvergence  for  general  meshes  with  local  re- 

finements 


Here  we  demonstrate  that,  for  all  practical  purposes,  the  17%-superconvergence 
regions  can  be  determined  for- any  element  in  the  interior  of  any  grid  with  local  refine¬ 
ments.  Let  us  consider  the  locally-refined  mesh  shown  in  Fig.  13a  and  let  us  assume 
that  we  are  interested  to  determine  the  17%-superconvergence  regions  for  the  deriva¬ 
tives  in  the  elements  shown  shaded  gray  in  the  interior  of  the  patch  of  interest  (which 
is  shown  with  thick  perigram  in  Fig.  13a).  In  order  to  find  the  17%-superconvergence 
regions,  the  patch,  which  is  shown  in  13b,  was  extracted  from  the  mesh  and  completed 
into  a  periodic  mesh  T  in  the  square-cell  c,  as  shown  in  Fig.  13c.  In  Fig.  14a  we 
show  the  regions  of  17%-superconvergence,  ,T)  for  17%  =  5%,  10%,  25%. 


In  Fig.  14b  we  show  the  bands  of  17%  relative  error  («;t,T^)  which  were  obtained 


by  solving  the  Neumarn  problem  in  the  mesh  shown  in  Fig.  13a  with  data  consistent 
with  the  harmonic  solution  v(zj,Z3)  =  sin(irxj)sinh(irx2).  It  can  be  observed  that  the 
regions  of  17%-superconvergence  are  practically  the  szmie  as  the  bands  of  17%  relative 
error.  (Here  we  observe  that  there  is  only  one  monomial  of  degree  (p-t- 1)  with  non-zero 
relative  error  in  the  Xj-derivative  in  the  elements  of  interest). 

We  also  repeated  the  above  computations  with  the  exact  solution  (Gaussian  hill 
function)  u(xi,X3)  =  with  fc  =  10  (note  that  this  solution  is  not 

harmonic).  The  center  of  the  Gaussian  hill  (a,  b)  is  allowed  to  take  values  of  (0.5, 
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0.5),  (-0.65,  -0.3)  and  (0.25,  0.35).  In  Figs.  15a,  16a  and  17a  we  show  the  bands 
of  i;%-relative  error  H^(it;r,Tj^)  obtained  by  solving  the  Neumann  boundary- value 

problem  in  the  mesh  Tf^  shown  in  Fig.  13a.  In  Figs.  15b,  16b  and  17b  we  give  the 
bands  of  17%-relative  error  ,7*)  obtained  by  solving  the  periodic  boundary- 


value  problem  over  the  periodic  mesh  T  shown  in  Fig.  13c  with  data  consistent  with  the 
Taylor  series  expansion  Q  of  the  Gaussian  hill  function.  It  can  be  observed  that  the  17%- 
relative  error  bands  obtained  from  the  actual  finite  element  solution  are,  for  all  practical 
purposes,  identical  to  those  obtained  from  the  periodic  boundary-value  problem  using 
the  Taylor  series  expansion.  In  Fig.  18  we  show  the  i/%-superconvergence  regions 
(G^;‘r,T)  for  17%  =  10%,  30%,  50%.  In  the  element  with  constrained  nodes  the 

17%-superconvergence  regions  for  the  class  of  general  solutions  and  17%  <  50%,  do  not 
include  the  center  of  the  element  which  is  the  Gauss- Legendre  point  for  p  =  1. 


5.6  77%-superconvergence  regions  for  a  complex  grid  of  skewed 
quadrilaterals  generated  by  a  mesh  generator 

We  will  now  give  an  example  of  finding  the  17%-superconvergence  regions  in  the 
interior  of  a  complex  grid  generated  by  commadal-like  mesh-generator.  We  considered 
the  mesh  of  skewed  quadrilaterals  shown  in  Fig.  19a.  We  were  interested  to  find  the 
17%-superconvergence  regions  in  the  elements  which  are  shown  shaded  in  Figs.  19a,  20a. 
In  each  case  we  considered  the  mesh-patches  which  include  three  mesh-layers  around 
the  elements  of  interest  as  shown  in  Figs.  19b,  20b.  These  mesh-patches  were  scaled 
and  translated  into  the  square-cell  c  :=  (—0.5, 0.5)  x  (—0.5, 0.5)  and  each  mesh-patch 
was  completed  into  a  periodic  mesh  (resp.  7^  )  as  shown  in  Fig.  19c  (resp.  Fig.  20c). 

In  order  to  demonstrate  that  the  results  for  the  17%-superconvergence  regions  ob¬ 
tained  from  the  methodology  are  also,  for  all  practical  purposes,  valid  in  the  actual 
mesh  we  performed  the  following  computations.  We  computed  the  17%-relative  error 
regions  with  the  error  obtained  by  solving  the  following  problems: 

PI:  Neumann  boundary- value  problem  (with  data  consistent  with  the  exact  solution 
u)  solved  using  the  domain  and  the  mesh  (actual  grid)  shown  in  Fig.  19a. 

P2:  Neumann  boundary- value  problem  (vrith  data  consistent  with  the  Taylor  series 
expansion  Q  of  the  exact  solution  u  about  the  center  of  the  mesh-patch  of  interest) 
solved  using  the  domain  and  the  mesh  (actual  grid)  shown  in  Fig.  19a. 

P3:  Periodic  boundary-value  problem  (with  data  consistent  with  the  Taylor  series 
expansion  Q  of  the  exact  solution  u  about  the  center  of  the  mesh- patch  of  interest) 
solved  using  the  periodic  dommn  and  mesh  (periodic  grid  obtained  by  freezing 
the  periodicity)  shown  in  Figs.  19c  and  20c. 
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We  computed  the  17%-relative  error  regions  for  the  follovring  solutions: 


(5.5o)  h  =  l,...,4 


(5.56) 


tv' 


W(z)  = 


o*  -  (2  +  4*^)’  o*  +  (z  +  4*^) 


(*) 


/  1  1  X 

P  Vo*  —  1  o’  + 


where  2  =  +  ias,  and  we  let  a  =  2  and  2^'^  =  |(1  + 1) ,  2©*^  =  |(— 1  +  0 »  4*^  = 

|(— 1  —  t)  ,  4*^  =  5(1  “0*  The  exact  solutions  given  above  are  analytic  in  the  domain 
and  are  t3rpical  solutions  of  elliptic  boundary- value  problems  in  two  dimensions  (see  also 
[25]).  We  computed  the  regions  of  25%  relative  error  in  the  shaded  elements  shown  in 
Figs.  19b,  20b  and  determined  the  25%-relative  error  region  which  is  common  to  all  four 
solutions  T,  T^^),  i  =  1, 2.  In  Figs.  21a,  21b,  21c  (resp.  21d,  21e,  21f)  we 
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show  the  25%-common  regions  (resp.  ^’#?({tvW}J_i;r,T;j)) 

in  the  elements  shown  in  Fig.  19a,  computed  by  solving  the  problems  PI,  P2,  P3, 
respectively.  It  can  be  observed  that  the  25%-common  regions,  computed  from  the 
three  problenu,  are  identical,  for  aU  practical  purposes.  In  Fig.  21g  (resp.  21h)  we 
show  the  17%-superconvergence  regions  (resp.  \Tyf\)).  Note 
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that  the  25%-superconvergence  regions  lie  entirely  within  the  25%-common  regions. 

Similarly  in  Figs.  22a,  22b,  22c  (resp.  22d,  22e,  22f)  we  show  the  25%-common 
regions  in  the  elements  shown 

in  Fig.  20a,  computed  by  solving  the  problems  PI,  P2,  P3,  respectively.  In  Fig.  22g 
(resp.  Fig.  22h)  we  show  the  17%-superconvergence  regions  r,  T^)  (resp. 

(Q^;  r,  7])).  Again  we  note  that  the  25%-superconvergence  regions  are  well  within  the 
corresponding  25%-common  regions. 


From  the  results  in  Sections  5.5  and  5.6  we  make  the  following  observations: 

(i)  The  17%-relative  error  regions  computed  &om  the  exact  solution  and  the  actual  or 
its  Taylor  series  expansion  using  the  periodic  grid  are,  for  all  practical  purposes, 
identical. 


(ii)  Since  we  have  studied  f7%-superconvergence  for  smooth  solutions  the  mesh  was 
allowed  to  be  quite  arbitrary. 

(iii)  The  locations  of  the  f7%-superconvergence  regions,  for  small  tf,  in  the  elements 
cannot  be  predicted  based  on  intuition.  In  many  instances  they  are  quite  far 
from  the  Gauss-Legendre  points  and  hence  the  Gauss-Legendre  points  may  not 
be  suitable  sampling  points. 
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5.7  The  pollution  effect  and  f7%*superconvergence  for  singu¬ 
lar  solutions 

In  this  Section  we  give  examples  which  show  the  effect  of  the  pollution  on  the 
superconvergence  of  finite  element  solutions.  In  particular  we  demonstrate  that: 

1.  In  finite  element  approximations  of  nngular  solutions  with  quasiuniform  meshes 
the  pollution  error  may  be  significant  (depending  on  the  strength  of  the  singu¬ 
larity,  the  polynomial  degree  of  the  elements  and  the  region  of  interest)  and  the 
i]i%-superconvergence  regions  may  not  exist,  asymptotically. 

2.  For  meshes  which  are  nearly  equilbrated  in  the  energy-norm  the  pollution-effect 
is  insignificant  and,  for  aU  practical  purposes,  we  can  have  17%-superconvergence 
of  the  derivatives  (as  predicted  by  the  methodology  of  Section  4)  in  any  interior 
element  of  the  mesh. 


To  show  the  pollution  effect  for  meshes  of  bilinear  elements  we  considered  the  L- 
shaped  domain  shown  in  Fig.  23  with  boundary  conditions  compatible  with  the  exact 
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solution  ^(r,^)  =  r“sin(a5).  For  a  =  -  (resp.  o  =  -)  homogeneous  Dirchlet 
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boundary-conditions  (resp.  homogeneous  mixed  boundary-conditions)  were  applied 
on  the  edges  of  the  reentrant  comer;  in  the  remaining  edges  of  the  boundary  Neu¬ 
mann  boundary-conditions  compatible  with  the  exact  solution  were  imposed.  Note 
that  the  theoretical  analysis  shows  (see  [24],  [26]  and  the  references  therein)  that  for 
quasiuniform  meshes  of  bilinear  elements  (p  =  1)  the  pollution  effect  is  significant 

(asymptotically)  for  a  <  -.  Thus  as  the  mesh  is  refined  uniformly  we  expect  to  see 
the  following: 

2 

(i)  For  a  =  -  the  q%-superconvetgence  regions  exist  asymptotically,  as  predicted 

O 

by  the  methodology. 

(ii)  For  a  =  ^  there  are  no  q%-superconvergence  regions,  asymptotically,  due  to  the 
pollution  effect. 


We  now  proceed  to  the  description  of  the  numerical  results. 

We  meshed  the  X-shaped  domain  with  uniform-mesh  of  bilinear  elements  and  con¬ 
sidered  the  mesh-patches  in  the  subdomains  shown  in  Fig.  24a  (subdomain  A)  and 
Fig.  24b  (subdomain  B).  In  Figs.  25a,  25b  and  25c  (resp.  Fig.  25d,  25e  and  25f)  we 
show  the  bands  (u;r,  T^)  in  subdomain  A  (resp.  subdomain  B)  for  the  uniform 

h  h  .  2 

meshes  with  element-size  h,  —  and  — ,  respectively,  for  a  =  -r .  It  can  be  observed  that 

«  4  o 
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the  17%-baiids  for  exist  asymptotically  and  coincide  with  the  bands  predicted  by 

the  methodology  of  Section  4.  In  Figs.  26a,  26b  and  26c  (resp.  Figs.  26d,  26e  and  26f) 

we  show  the  bands  H^(it;r,T^)  for  a  —  it  can  be  observed  that  the  i7%-bands 

disappear  as  the  mesh-size  h  tends  to  zero.  Tht  reason  for  the  disappearance  of  the 
bands  is  that  the  error  is  uniformly  large  in  the  entire  element  (see  also  [26]). 

We  now  consider  the  class  of  grids  which  are  nearly  equilibrated  in  the  energy-norm 
which  are  obtained  by  refining  adaptively  a  coarse  mesh  of  squares.  In  Fig.  27a  (resp. 

a  2 

Fig.  28a)  we  show  a  nearly  equilibrated  grid  for  the  exact  solution  u(r,  9)  =  r^  sin(-d) 

O 

and  mark  the  elements  of  interest  in  Patch  P  (resp.  Patch  Q).  The  periodic  meshes 
obtained  by  completing  these  patches  are  shown  in  Fig.  27b  and  Fig.  28b  for  patches  P 
and  Q,  respectively.  In  Fig.  29a  and  (resp.  Fig.  30a)  we  show  the  bands  (u;r,T^) 

obtained  from  the  solution  of  the  finite  element  problem  on  the  actual  mesh  and  in 
Fig.  29b  (resp.  Fig.  30b)  we  show  the  corresponding  q%-bands  obtained  from  the 
periodic  boundary- value  problem  using  Taylor  series  expansion  on  the  mesh  shown  in 
Fig.  27b  (resp.  Fig  28b),  respectively. 

We  also  considered  the  nearly  equilibrated  grid  for  the  exact  solution  u(r,0)  = 

rT  sin(-d)  with  appropriate  boundary  conditions,  as  discussed  earlier.  The  mesh  and 
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the  elements  of  interest  are  shown  in  Fig.  31a  and  Fig.  32a.  The  corresponding  periodic 
meshes  for  these  patches  are  shown  in  Fig.  31b  and  Fig.  32b.  We  computed  the  bands 
(u;  r,  Tf^)  using  the  finite  element  solution  in  the  actual  mesh  shown  in  Fig.  31a  (or 

Fig.  32a)  and  the  bands  5*'^  (Q;  r,  T)  from  the  finite-element  solution  of  the  periodic 

boundary-value  problem  using  the  mesh  shown  in  Fig.  31b  (or  Fig.  32b).  These  are 
shown  in  Fig.  33a  and  Fig.  33b  (resp.  Fig.  34a  and  Fig.  34b)  for  the  patch  shown  in 
Fig.  31a  (resp.  Fig.  32a).  It  can  be  seen  that  the  corresponding  bands  are  identical. 

Hence  we  conclude  that  when  the  mesh  is  adaptive  (nearly  equilibrated  in  the  energy 
norm)  the  17%-superconvergence  regions,  which  are  determined  using  the  methodology 
of  Section  4,  are  valid  for  any  interior-patch  of  the  mesh. 


6  Summary  of  conclusions 

We  presented  the  new  concept  of  17%-superconvergence  which  enables  us  to  deter¬ 
mine  the  points  of  least  error  for  the  derivatives  in  elements  in  the  interior  of  general 
meshes  of  quadrilaterals.  We  conducted  a  numerical  study  to  determine  the  applica¬ 
bility  of  the  new  concept  in  practical  grids  and  we  arrived  at  the  following  conclusions: 

1.  The  regions  of  T7%-superconvergence  in  any  mesh  can  be  determined  by  using  a 
computer-based  approach.  This  approach  takes  directly  into  account  the  topol- 
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ogy  of  the  mesh,  the  degree  p  of  the  elements  ud  the  nature  of  the  solution 
(harmonic  or  general). 

2.  The  classical  superconvergence  points  (i.e.  the  0%-superconvergence  points)  in 
an  element  are  very  sensitive  to  the  local  geometry  of  the  mesh  and  in  some  cases 
they  may  not  exist. 

3.  The  17%-superconvergence  regions  for  the  derivatives  in  an  interior  element  are 
robust  with  respect  to  the  topology  of  the  mesh  outside  the  immediate  neigh¬ 
borhood  of  the  element  and  can  always  be  determined  for  any  class  of  smooth 
solutions  using  the  method  of  freezing  the  periodicity. 

4.  For  singular  solutions  computed  using  quasi-uniform  meshes,  the  pollution  error 
may  be  significant  and  q%-superconvergence  regions  may  not  exist  asymptotically 
for  these  meshes.  However,  for  meshes  which  are  nearly  equilibrated  in  the  energy 
norm,  the  pollution  error  in  the  derivatives  is  negligible  and  q%-superconvergence 
regions  for  the  derivatives  cui  be  determined  in  any  interior  element  of  the  mesh. 

5.  The  position  of  the  superconvergence  points  for  the  typical  local  topologies,  as 
reported  here,  can  be  used  if  the  mesh  generator  creates  meshes  with  a  small 
number  of  local  mesh-topologies  as  for  example  in  the  codes  using  constrained 
square  meshes  with  refinements  (or  mapped  meshes). 

The  methodology  presented  in  this  paper  can  be  also  employed  to  study  the  rf%- 
superconvergence  of  flux-quantities  in  finite  dement  approximations  of  elliptic  systems 
(f7%-superconvergence  of  the  gradient  of  the  displacement,  strain,  stress  in  elasticity 
etc.). 
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List  of  Figures 

Fig.  1.  ^-refinement  for  bilinear  and  biquadratic  quadrilaterals:  1-irregular  rule.  The 
active  nodes  are  shown  with  a  soUd  circle  and  the  constrained  nodes  are  indicated  by 
a  circle,  (a)  A  mesh-patch  of  bilinear  elements  showing  the  locations  of  active  and 
constrained  nodes;  (b)  A  mesh-patch  of  biquadratic  elements  showing  the  location  of 
active  and  constrained  nodes;  (c)  The  mesh-patch  of  bilinear  quadrilaterals  shown  in 
(a)  with  transition  elements;  (d)  The  mesh-patch  of  biquadratic  elements  shown  in  (b) 
with  transition  elements. 


Fig.  2.  Locally  periodic  grid:  (a)  An  X-shaped  domain  with  a  periodic-mesh  subdo¬ 
main  XT)  covered  exactly  by  a  periodic  array  of  h-cells.  The  perimeter  of  the 

periodic-mesh  subdomain  5(e°,  H)  is  shown  with  a  thick  line;  (b)  The  master-mesh  T 
in  the  periodic-mesh  subdomain  5(8!*’,  H°). 

Fig.  3.  Method  of  freezing  the  periodicity:  (a)  A  finite  element  mesh  of  skewed 
quadrilaterals  with  the  interior  mesh-patch  o/g  shown  shaded  gray  and  the  mesh-patch 
ufg  (which  includes  Ug  in  its  interior)  shown  with  thick-black  perigram;  (b)  Extraction 
of  the  mesh-patch  Ug  and  completion  to  a  periodic  mesh  in  a  square-cell;  (c)  The  finite 
element  mesh  assumed  by  the  method  of  freezing  the  periodicity. 

Fig.  4.  Typical  mesh  topologies  in  locally  refined  meshes:  (a)  Mesh-interface  with 
constrained  nodes;  (b)  Mesh-interface  with  transition  elements;  (c)  A  periodic  mesh 
with  two  levels  of  refinement. 


Fig.  5.  Mesh-interface  with  constrained  nodes:  Superconvergence  points,  (a)  Bilinear 
elements:  Zero-lines  of  ’  (thick-lines)  and  u---  ■  (thin-lines)  intersecting  at  the 


dxi 


points  of  simultaneous  superconvergence  for  the  x^-  and  Xj-derivatives  for  the  class  of 
harmonic  solutions;  (b)  Bilinear  elements:  Zero-lines  of  ^  ’  (thick-lines)  and  ^ 

C/X2 

(thin-lines)  intersecting  at  the  points  of  superconvergence  for  the  X2-derivative  for  the 
class  of  general  solutions;  (c)  Biquadratic  elements:  Zero-lines  of  —  (thick-lines) 

CrXj 

and  -^-2-  (thin-lines)  and  their  intersections.  Note  that  does  not  vanish  (in  the 

elements  of  size  —  but  is  very  small)  at  these  points  which  are  0.001%  superconvergence 
2 

points. 


Fig.  6.  Rate  of  convergence  at  the  superconvergence  points:  Sequence  of  meshes  with 
a  mesh-interface.  Meshes  with  constrained  nodes  along  the  interface.  The  cell-size  in 
the  meshes  is  (a)  h  =■  0.25;  (b)  h  =  0.125;  (c)  h  =  0.0625;  (d)  h  =  0.03125. 

Fig.  7.  Rate  of  convergence  at  the  superconvergence  points:  Sequence  of  meshes  with 
a  mesh-interface.  Meshes  with  transition-elements  along  the  interface.  The  cell  size  in 
the  meshes  is  (a)  h  =  0.25;  (b)  h  =  0.125;  (c)  h  =  0.0625;  (d)  h  =  0.03125. 
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Fig.  8.  Relative  error  at  the  superconvergence  points:  An  adaptive  mesh  of  bilinear 
elements  generated  for  the  exact  solution  u(xj,X3)  =  used  for  the 

study,  (a)  An  interior  patch  of  elements  with  a  mesh-interface  (Patch  I)  in  the  finite 
element  mesh;  (b)  The  mesh-patch  I  with  the  elements  of  interest  enumerated;  (c)  An 
interior  patch  of  elements  with  a  small  region  with  a  mesh-interface  (Patch  II);  (d) 
The  mesh-patch  II  with  the  elements  of  interest  enumerated. 

Fig.  9.  17%-superconvergence  in  periodic  meshes  with  local  refinements:  Contours 
of  0%-error  and  17%-superconvergence  regions  for  the  harmonic  monomials  of  degree 
(p-l- 1).  The  elements  shown  shaded  gray  in  the  mesh  (to  the  right)  are  considered,  (a) 
C^{Q^\r,f)  and  1Vt{Q^)T,f)  for  p  =  1;  (b)  C^(Q®;T,r)  and 

for‘p  =  2;  (c)  (C^r.f)  and  ((2^;r,f)  for  p  =  1;  (d)  C®*  ((2^;r,f)  and 

(C^r,f)  for  p  =  2. 

Fig.  10.  f/%-superconvergence  in  periodic  meshes  with  local  refinements:  Mesh  with 
transition  elements  at  the  mesh  interface.  Contours  of  0%-error  and  17%-superconvergence 
for  the  class  of  harmonic  monomials  of  degree  (p-l-1).  (a)  t,  f)  and  TJ®*  (G^;  t,  T) 

for  p  =  1;  (b)  C®*  (G^;r,r)  and  TVt  {Q^]T,f)  for  p  =  1;  (c)  C®5  (G®;r,r)  and 

9m2  Ai}  3mi 

'RP’SL  (G^JTjT)  for  p  =  2;  The  thin  line  is  the  0%-error  contour  for  the  monomial  Qf 

and  the  thick  line  is  the  0%-error  contour  for  the  monomial  Qf .  Note  that  in  the 
triangular  subelements  we  have  superconvergence  points  and  in  the  square  elements 
we  have  bands  of  superconvergence. 

Fig.  11.  f7%-superconvergence  in  periodic  meshes  with  local  refinements:  Mesh  with 
transition  elements  at  the  mesh  interface.  Contours  of  0%  error  and  q%-superconvergence 
regions  for  the  class  of  general  monomials  of  degree  (p  -|- 1). 

(a)  C®’^  (G^;r,r)  and  TV'l  (Q®;r,r)  for  p  =  1; 

(b)  C®^  (G^;r,r)  and  TVl  {Q^]r,f)  for  p  =  1; 

992 

(c)  C°^  (Q^;  r,  f)  and  TV'l  (G^;  r,  f )  for  p  =  2; 

(d)  C®^  (G^;r,r)  and  TV'l  for  p  =  2. 

Note  that  for  p  =  2,  there  are  no  regions  of  q%-superconvergence  for  17%  <  5%. 

Fig.  12.  f7%-superconvergence  in  periodic  meshes  with  local  refinements  as  shown 
in  Fig.  4c.  Contours  of  0%  error  and  77%-superconvergence  r^ons  for  the  class  of 
harmonic  monomials  of  degree  (p-M).  (a)  C®^  and  (G^;‘r,T)  forp  =  1; 

(b)  (G";  r,  T)  and  (G^;  r,  f)  for  p  =  2;  (c)  C°’‘  (G^;  r,  T)  and  7V>^  (G^;  r,  f) 

SSj  d«2 

for  p  =  1;  (d)  C®(S  (G^;t, T)  and  {Q^]t,T)  for  p  =  2.  The  thin  line  corresponds 

092  0«2 

to  the  0%-error  contour  for  the  monomial  Qf  and  the  thick  line  corresponds  to  the 
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0%-contoiii  for  the  monomial  Qf.  Note  that  in  some  elmnents,  the  17%  at  the  Ganss- 
Legendre  point  is  more  than  25%. 

Fig.  13.  f7%'8uperconvergence  for  general  meshes  with  local  refinements:  Adap¬ 
tive  mesh  of  bilinear  square  elements  generated  using  the  exact  solution  tt(zi,X3)  = 
g-ii((»i-a)»+(»a-fc)»]^  jg  _  finite  dement  mesh  Tjj  in  the  domain  ft.  The 

interior  mesh-patch  of  interest  is  shown  with  a  thick  perigram;  (b)  The  interior  mesh- 
patch  of  interest  (shown  in  (a))  is  shown  with  the  elements  in  which  superconvergence 
is  studied  (shown  with  gray  shading);  (c)  The  periodic  mesh  f  obtained  by  completing 
the  patch  of  interest.  The  patch  of  elements  in  which  superconvergence  is  studied  are 
shown  with  gray  shading. 

Fig.  14.  q%-superconvergence  for  general  meshes  with  local  refinements:  Adaptive 
mesh  of  bilinear  square  elements  shown  in  Fig.  13a.  (a)  q%-superconvergence  regions 
^’'5  (G®;  Ti,T)  for  the  harmonic  solutions  in  the  elements  marked  gray  in  the  peri- 

odic  mesh  T  shown  in  Fig.  13c;  (b)  Bands  of  17%-relative  error  for  the  exact  solution 
u(zj,X3)  =  sin(TXj)sinh(wx2)  in  the  elements  marked  gray  in  the  finite  element  mesh 
Tf^  as  shown  in  Fig.  13b.  Note  that  the  17%-superconvergence  regions  shown  in  (a)  are 
very  similar  to  the  bands  of  17%-relative  error  shown  in  (b). 

Fig.  15.  17%-superconvergence  for  general  meshes  with  local  refinements:  Adap¬ 
tive  mesh  of  bilinear  elements  shown  in  Fig.  13a.  (a)  Bands  of  17%'relative  error 
for  the  exact  solution  where  k  =  10  and 

(a,b)  =  (0.5,  0.5)  computed  by  solving  the  Neumann  problem  in  the  mesh  shown 
in  Fig.  13a.  The  bands  of  q%-relative  error  shown  for  the  elements  marked  gray  in 
Fig.  13b;  (b)  Bands  of  »/%-relative  error  (Q;  r,  f)  for  the  Taylor  series  expansion 

of  u  computed  by  solving  the  periodic  boundary-value  problem  in  the  periodic  mesh 
T  shown  in  Fig.  13c.  The  bands  of  t7%-relative  error  shown  for  the  elements  marked 
gray  in  Fig.  13c.  Note  that  the  bands  of  17%-relative  error  computed  from  the  finite 
element  solution  are  very  similar  to  the  bands  of  17%-relative  error  computed  from  the 
periodic  boundary-value  problem. 

Fig.  16.  17%-superconvergence  for  general  meshes  with  local  refinements:  Adaptive 
mesh  of  bilinear  square  elements  shown  in  Fig.  13a:  (a)  Bands  of  17%-relative  error 
{u]T,Tf^)  for  the  exact  solution  ^(xj,®,)  =  where  fc  =  10  and 

(a,  6)  =  (-0.65,  -0.30)  computed  by  solving  the  Neumann  problem  in  the  mesh  shown 
in  Fig.  13a.  The  bands  of  17%-relative  error  shown  for  the  elements  marked  gray  in 
Fig.  13b;  (b)  Bands  of  17%-relative  error  (Qi^iT)  for  the  Taylor  series  expansion 

of  u  computed  by  solving  the  periodic  boundary-value  problem  in  the  periodic  mesh 
f  shown  in  Fig.  13c.  Note  that  the  bands  obtained  from  the  above  figures  are  very 
similar. 
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Fig.  17.  )7%-siiperconTergeiice  for  general  meshes  with  local  refinements:  Adaptive 
mesh  of  bilinear  square  elements  shown  in  Fig.  13a:  (a)  Bands  of  17%-relative  error 

B”*  for  the  exact  solution  where  ib  =  10  and 

wrr 

(a,  b)  =  (0.25,  0.35)  computed  by  solving  the  Neumann  problem  in  the  mesh  shown 
in  Fig.  13a.  The  bands  of  q%-relative  error  shown  for  the  elements  marked  gray  in 
Fig.  13b;  (b)  Bands  of  17%-relative  error  (Q;  r,  T)  for  the  Taylor  series  expansion 

of  u  computed  by  solving  the  periodic  boundary-value  problem  in  the  periodic  mesh 
T  shown  in  Fig.  13c.  Note  that  the  bands  of  17%-relative  error  given  in  the  above  two 
figures  are  very  similar. 

Fig.  18.  17%-superconvergence  for  general  meshes  with  local  refinements:  Adaptive 
mesh  of  bilinear  square  elements  shown  in  Fig.  13b.  Regions  of  17%-superconvergence 
(Q^;r,T)  for  the  general  solutions  in  the  elements  marked  gray  in  the  periodic 

mesh  T  (shown  in  Fig.  13c).  Note  that  regions  are  shown  for  17%  =  10%,  30%,  50% 
and  that  in  some  elements  the  regions  do  not  exist  for  17%  <  25%. 

Fig.  19.  17%-superconvergence  for  a  grid  generated  by  a  mesh  generator:  Mesh  of 
skewed  quadrilaterals  generated  by  converting  a  mesh  of  triangles,  (a)  The  finite 
element  mesh  in  the  domain.  The  patch  of  elements  in  which  superconvergence  is 
studied  are  marked  with  gray  shading.  The  patch  surrounding  these  elements  is  shown 
with  a  thick  perigram.  (b)  The  patch  of  elements  shown  in  (a);  (c)  The  periodic  mesh 
T  generated  by  completing  the  patch  shown  in  (b). 

Fig.  20.  17%-superconvergence  for  a  grid  generated  by  a  mesh  generator:  Mesh  of 
skewed  quadrilaterals  generated  by  converting  a  mesh  of  triangles,  (a)  The  finite 
element  mesh  in  the  domain.  The  patch  of  elements  in  which  superconvergence  is 
studied  are  shown  with  gray  shading.  The  mesh  patch  surrounding  these  elements  is 
shown  with  a  thick  perigram;  (b)  The  patch  of  elements  marked  in  (a);  (c)  The  periodic 
mesh  T  generated  by  completing  the  patch  shown  in  (b). 

Fig.  21.  17%-superconvergence  for  a  grid  of  skewed  quadrilaterals  generated  by  a  mesh 
generator:  Elements  shown  with  gray  shading  in  Fig.  19a.  Contours  of  25%-relative 
error  for  the  four  solutions  given  in  Section  5.6.  The  regions  of  25%-relative  error 
which  are  common  to  the  four  solutions  are  shaded  gray.  Bilinear  elements  (p  = 
1).  (a)  Contours  of  25%-relative  error  and  25%-common  regions  for  the  ®j-derivative 
computed  by  solving  the  problem  PI  (Neumann  problem  with  data  consistent  with 
the  exact  solution);  (b)  Contours  of  25%-relative  error  and  25%-common  regions  for 
the  ®j-derivative  computed  by  solving  the  problem  P2  (Neumann  problem  with  data 
consistent  with  Taylor  series  expansion);  (c)  Contours  of  25%-relative  error  and  25%- 
common  regions  for  the  x^-derivative  computed  by  solving  the  problem  P3  (Periodic 
boundary-value  problem  vnth  Taylor  series  expansion);  (d)  Contours  of  25%-relative 
error  and  25%-common  regions  for  the  ®j-derivative  computed  by  solving  problem 
PI  (Neumann  problem  over  the  entire  domain);  (e)  Contours  of  25%-relative  error 
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&nd  25%-common  regions  for  the  Xj-derivative  computed  by  solving  the  problem  P2 
(Neumann  problem  with  data  consistent  with  Taylor  series  expansion);  (f)  Contours  of 
25%-relative  error  and  25%-common  regions  for  the  Xj-derivative  computed  by  solving 
the  problem  P3  (Periodic  boundary-value  problem  with  Taylor  seriM  expansion);  (g) 
17%-superconvergence  regions  for  the  Xj-derivative  for  the  harmonic 

solutions  of  degree  2.  (h)  i/%-superconvergence  regions  {.Q,^\r,f)  for  the  x,- 

derivative  of  the  error  for  the  harmonic  solutions  of  degree  2.  Note  that  Figures  (a), 
(b)  and  (c)  (resp.  Figures  (d),  (e)  and  (f))  are,  for  all  practical  purposes,  identical. 

Fig.  22.  17%-superconvergence  for  a  grid  of  skewed  quadrialterals  generated  by  a  mesh 
generator:  Elements  shown  with  gray  shading  in  Fig.  20a.  Contours  of  25%-relative 
error  for  the  four  solutions  given  in  Section  5.6.  The  regions  of  25%-relative  error 
which  are  common  to  the  four  solutions  are  shaded  gray.  Bilinear  elements  (p  = 
1).  (a)  Contours  of  25%-relative  error  and  25%-common  regions  for  the  x^-derivative 
computed  by  solving  the  problem  PI  (Neumann  problem  with  data  consistent  with 
the  exact  solution);  (b)  Contours  of  25%-relative  error  and  25%-common  regions  for 
the  Xj -derivative  computed  by  solving  the  problem  P2  (Neumann  problem  with  data 
consistent  with  Taylor  series  expansion);  (c)  Contours  of  25%-relative  error  and  25%- 
common  regions  for  the  x^-derivative  computed  by  solving  the  problem  P3  (Periodic 
boimdary-value  problem  with  Taylor  series  expansion);  (d)  Contours  of  25%-relative 
error  and  25%-common  regions  for  the  x,-derivative  computed  by  solving  problem 
Pi  (Neumann  problem  over  the  entire  domain);  (e)  Contours  of  25%-relative  error 
and  25%-common  regions  for  the  x, -derivative  computed  by  solving  the  problem  P2 
(Neumann  problem  with  data  consistent  with  Taylor  series  expansion);  (f)  Contours  of 
25%-relative  error  and  25%-common  regions  for  the  Xj-derivative  computed  by  solving 
the  problem  P3  (Periodic  boundary-value  problem  with  Taylor  series  expansion);  (g) 
17%-superconvergence  regions  (Q^;r,T)  for  the  Xj-derivative  for  the  harmonic 

solutions  of  degree  2.  (h)  q%-superconvergence  regions  (Q^;t,T)  for  the  Xj- 

derivative  of  the  error  for  the  harmonic  solutions  of  degree  2.  Note  that  Figures  (a), 
(b)  and  (c)  (resp.  Figures  (d),  (e)  and  (f))  are,  for  all  practical  purposes,  identical. 

Fig.  23.  Pollution  effect  and  q%-superconvergence  for  singular  solutions:  X-shaped 
domain. 

Fig.  24.  Pollution  effect  and  77%-superconvergence  for  singular  solutions:  X-shaped 
domain  meshed  with  a  uniform  mesh  of  square  elements,  (a)  Subdomain  A;  (b)  Sub- 
domain  B. 


Fig.  25.  Pollution  effect  and  q%-superconvergence  for  singular  solutions:  X-shaped 

3  .  2 

domain  meshed  with  a  uniform  mesh  of  squares;  exact  solution  u(r,0)  =  r?  sin -d; 

O 

bilinear  elements  (p  =  1).  Bands  of  17%- relative  error  for  the  Xj-derivative  (u;  r,  T^) 
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for  the  elements  in:  (a)  Subdomain  A  with  element  size  h;  (b)  Subdomain  A  with 

h  .  h  . 

element  size  — ;  (c)  Subdomain  A  with  element  size  — ;  (d)  Subdomain  B  with  element 
2  ^4  ^ 

size  h]  (e)  Subdomain  B  with  element  size  — ;  (f)  Subdomain  B  with  element  size  — . 

A  4 

Note  that  as  the  element  size  is  decreased  the  bands  converge  to  the  bands  predicted 
by  the  methodology  of  Section  4. 

Fig.  26.  Pollution  effect  and  i/%-superconTergence  for  singular  solutions:  Z-shaped 

domain  meshed  with  a  uniform  mesh  of  squares.  Exact  solution  u(r,  6)  =  ri  sin  -9, 

o 

p  =  1.  Bainds  of  17%-relative  error  for  the  acj-derivative  for  the  elements 

h 

in:  (a)  Subdomain  A  with  element  size  h;  (b)  Subdomain  A  with  element  size  — ; 

A 

(c)  Subdomain  A  with  element  size  — ;  (d)  Subdomain  B  with  element  size  h;  (e) 

h  h 

Subdomain  B  with  element  size  — ;  (f)  Subdomain  B  with  element  size  — .  Note  that 

A  4 

as  the  element  size  is  decreased  the  bands  of  25%-relative  error  disappear. 

Fig.  27.  Pollution  effect  and  17%-superconvergence  for  singular  solutions:  Z-shaped 

2  2 

domain;  exact  solution  =  risin-^;  bilinear  elements  (p  =  1).  (a)  Adaptive 

mesh  of  square  elements  (nearly  equilibrated  in  the  energy  norm).  The  patch  of  ele¬ 
ments  in  which  superconvergence  is  studied  are  shown  with  gray  shading.  The  patch 
surrounding  these  elements  is  marked  with  a  thick  perigram  (Patch  P);  (b)  Periodic 
mesh  T  obtained  by  completing  the  interior  patch  shown  in  (a). 

Fig.  28.  Pollution  effect  and  f7%-superconvergence  for  singular  solutions:  Z-shaped 

2  2 

domain;  exact  solution  v,(r,6)  =  r^sin— d;  bilinear  elements  (p  =  1).  (a)  Adaptive 

o 

mesh  of  square  elements  (nearly  equilibrated  in  the  energy  norm).  The  patch  of  ele- 
ments  in  which  superconvergence  is  studied  are  shown  with  gray  shading.  The  patch 
surrounding  these  elements  is  marked  with  a  thick  perigram  (Patch  Q);  (b)  Periodic 
mesh  T  obtained  by  completing  the  interior  patch  shown  in  (a). 

Fig.  29.  Pollution  effect  and  77%-superconvergence  for  singular  solutions:  Z-shaped 

2  2 

domain;  exact  solution  v,{r,6)  =  rr  sin bilinear  elements  (p  =  1).  The  patch  of 

o 

elements  shown  with  gray  shading  in  Fig.  27a  (Patch  P).  (a)  Bands  of  q%-relative  error 
for  the  ®i-derivative  {u•,T,Tf^)  computed  by  solving  the  finite  element  problem  for 

elements  in  Patch  P;  (b)  Bands  of  77%-relative  error  for  the  a^-derivative  {Q\  r,  T) 

computed  by  solving  the  periodic  boundary-value  problem  with  data  consistent  with 
the  Taylor  series  expansion  of  u. 
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Fig.  SO.  Pollution  effect  and  i7%-8uperconTergence  for  singular  solutions:  i^shaped 

3  2 

domain;  exact  solution  v(r,d)  =  rT sin-0;  bilinear  elements  (p  =  1).  The  patch  of 
elements  shown  with  gray  shading  in  Fig.  28a  (Patch  Q).  (a)  Bands  of  17%-relatiTe  error 
for  the  z^-derivative  (u; r^T^)  computed  by  solving  the  finite  element  problem  for 

elements  in  Patch  Q;  (b)  Bands  of  17%-relative  error  for  the  x^-derivative  B^(Q;  r,  T) 

computed  by  solving  the  periodic  boundary-value  problem  with  data  consistent  with 
the  Taylor  series  expansion  of  u. 

Fig.  31.  Pollution  effect  and  i7%-superconvergence  for  singular  solutions:  X-shaped 

domain;  exact  solution  =  r^sin— 0;  bilinear  elements  (p  =  1).  (a)  Adaptive 

mesh  of  square  elements  (nearly  equilibrated  in  the  energy  norm).  The  patch  of  ele¬ 
ments  in  which  superconvergence  is  studied  are  shown  with  gray  shading.  The  patch 
surrounding  these  elements  is  marked  with  a  thick  perigram.  (b)  Periodic  mesh  T 
obtained  by  completing  the  interior  patch  shown  in  (a). 

Fig.  32.  Pollution  effect  and  q%-superconvergencefor  singular  solutions:  X-shaped  do¬ 
main.  Exact  solution  it(r,  0)  =  ri  sin  -0.  p  =  1.  (a)  Adaptive  mesh  of  square  elements 

(nearly  equilibrated  in  the  energy  norm).  The  elements  in  which  superconvergence  is 
studied  are  shown  with  gray  shading.  The  patch  surrounding  these  elements  is  marked 
with  a  thick  perigram.  (b)  Periodic  mesh  T  obtained  by  completing  the  interior  patch 
shown  in  (a). 

Fig.  33.  Pollution  effect  and  i;%-superconvergence  for  singular  solutions:  X-shaped 
domain;  exact  solution  u(r,0)  =  r?  sin— 0;  bilinear  elements  (p  =  1).  The  patch  of 
elements  shown  with  gray  shading  in  Fig.  31a.  (a)  Bands  of  q%-relative  error  for  the 
Xj-derivative  (u;  r,  T^)  computed  by  solving  the  finite  element  problem;  (b)  Bands 

of  q%-relative  error  for  the  Xj -derivative  B’’^  (u;  r,  T)  computed  by  solving  the  periodic 

1 

boundary-value  problem  with  data  consistent  with  Taylor  series  expansion  of  ri  sin  -0. 

U 

Fig.  34.  Pollution  effect  and  q%-superconvergence  for  singular  solutions:  X-shaped 
domain;  exact  solution  u(r,  0)  =  r?  sin  -0;  bilinear  elements  (p  =  1).  Elements  shown 
with  gray  shading  in  Fig.  32a.  (a)  Bands  of  17%-relative  error  for  the  x^-derivative 
B’52  (“u;  “T,  T^)  computed  by  solving  the  finite  element  problem;  (b)  Bands  of  17%-relative 

error  for  the  Xi-derivative  B’S  (u;  r,  T)  computed  by  solving  the  periodic  boundary- 

*•1  ^ 

value  problem  with  data  consistent  with  Taylor  series  expansion  of  r?  sin  -0. 

o 
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Simultaneous  superconvergence  points  for  the  x^-  and  Xj-derivatives 
in  the  elements  near  the  mesh-interface  of  Fig.  4a 

Element 

p=  1 

p  =  2 

(  0.3449514150,  0.189895754) 

(-.2378069181,  -.5773502692) 
(-.5421972529,  .5773502692) 

(  .6071187677,  .5773502692) 

^2 

(-0.3449514150, 0.189895754) 

(  .2378069181, -.5773502692) 
(-.6071187677,  -.5773502692) 

(  .5421972529,  .5773502692) 

^3 

(  0.0000001171,  0.000000000) 

(-.6111972936,  -.5773502692) 

(  .5449056851,  -.5773502692) 

(  .5773502692,  .5773502692) 
(-.5773502692,  .5773502692) 

^4 

(-0.0000001171,  0.000000000) 

(  .6111972936,  -.5773502692) 
(-.5449056851,  -.5773502692) 

(  .5773502692,  .5773502692) 
(-.5773502692,  .5773502692) 

Table  1.  Superconvergence  points  in  the  elements  near  the  mesh-interface  with  con¬ 
strained  nodes:  Master  element  coordinates  of  the  points  for  the  elements  shown  in 
Fig.  4a.  The  points  of  simultaneous  superconvergence  for  the  Xj-  and  Xj-derivatives 
are  given.  Note  that  for  p  =  1,  the  superconvergence  points  in  Tj  and  Tj  are  not  close 
to  the  center  of  the  element  and  that  for  p  =  2,  there  are  only  three  superconvergence 
points  in  the  elements  Tj  and  Tj  .  The  points  for  p  =  2  are  T7%-superconvergence  points 
with  »/%  <  0.001%. 
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Superconvergence  points  for  the  Zi-,  Zj-derivative  in  the  elements 
near  the  mesh-interface  of  Fig.  4b;  linear  elements 


Element 

^(Q°;r,,f) 

W(C®;r.,r) 

(0.5000000000, 

0.5000000000) 

YES 

(0.2219334168, 

0.2291334168) 

NO 

(0.0000000000, 

0.6500546310) 

NO 

(0.7743219764, 

0.0000000000) 

NO 

^3 

(0.6500546310, 

0.0000000000) 

NO 

(0.0000000000, 

0.7743219764) 

NO 

(0.0000000000, 

0.9038367177) 

YES 

(-0.4759591794, 

0.0000004630) 

YES 

T-s 

(0.0000000000, 

0.9038367177) 

YES 

(0.4759591794, 

0.0000004630) 

YES 

(0.0000000000, 

0.8164965809) 

YES 

(0.0000000100, 

0.000000000) 

YES 

^7 

(0.0000000000, 

0.8164965809) 

YES 

(-0.0000000100, 

0.000000000) 

YES 

Table  2a.  Superconvergence  points  in  the  elements  near  the  mesh-interface  with 
transition  elements:  Master  element  coordinates  of  the  points  for  the  elements  shown  in 
Fig.  4b;  linear  elements.  The  superconvergence  points  for  the  triangular  sub-elements 
are  given  in  the  m2tster-element  coordinates  of  each  triangle.  The  points  which  are  also 
superconvergence  points  for  the  general  solutions  are  indicated  by  “yes”. 
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Superconvergence  points  for  the  x^-,  x^-derivative  in  the  elements 
near  the  mesh-interface  of  Fig.  4b;  quadratic  elements 

Element 

A^(G®;r.,f) 

T-i 

(0.6346980477, 

NO 

(0.3376461192, 

NO 

0.2686462265); 

0.3623896091); 

(0.2686462265, 

NO 

(0.0555703987, 

NO 

0.6346980477); 

0.0555703987); 

(0.3978751334, 

NO 

0.3978751334); 

(0.2960013625, 

NO 

0.2960013625); 

(0.2728756106, 

NO 

(0.0000000000, 

NO 

0.2305710433); 

0.1358925622); 

(0.6495869087, 

NO 

(0.0000000000, 

NO 

0.3168545738); 

0.7973549868); 

(0.2899495590, 

NO 

0.4545768446); 

(0.4561334539, 

NO 

0.2890038707); 

^3 

(0.2305710433, 

NO 

(0.1358925622, 

NO 

0.2728756106); 

0.0000000000); 

(0.3168545738, 

NO 

(0.7973549868, 

NO 

0.6495869087); 

0.0000000000); 

(0.4545768446, 

NO 

0.2899495590); 

(0.2890038707, 

NO 

0.4561334539); 

(0.5779800281, 

NO 

(1.0000000000, 

NO 

-0.1233185871); 

-0.7250525375); 

(-0.5783530068, 

NO 

(1.0000000000, 

NO 

0.0846847343); 

0.5795077940); 

(-0.5775012409, 

NO 

(-1.000000000, 

NO 

0.9724859376); 

0.5546779289); 

(-1.000000000, 

NO 

0.4051347137) 

Table  2b.  Superconvergence  points  for  the  elements  near  the  mesh-interface  with 
transition  elements:  Master  element  coordinates  of  the  points  for  the  elements  shown 
in  Fig.  4b;  quadratic  elements.  The  superconvergence  points  for  the  triangular  sub¬ 
elements  are  given  in  the  master-element  coordinates  of  each  triangle.  The  points  which 
are  also  superconvergence  points  for  the  general  solutions  are  indicated  by  “yes”.  Note 
that  there  are  no  superconvergence  points  for  the  class  of  general  solutions. 


Simultaneous  superconvergence  points  for  the 

Xi>  and  Xj-derivatives  for  the  mesh  of  Fig.  4c 

Bilinear  elements 

£lement 

^6 

(  1.0000000000, 1.0000000000) 

^8 

(  0.6000010117,  0.2000010310) 

(  0.2000001110,  0.6000000313) 

Tio 

(-1.0000000000, 1.0000000000) 

(-0.6000010117,-0.1999987998) 

(  0.6000010117,-0.1999987998) 

^13 

(  1.0000000000, 1.0000000000) 

(-0.2000001110,-0.6000000313) 

ns 

(-1 .0000000000,- 1 .0000000000) 

ne 

(-1.0000000000, 1.0000000000) 

nr 

(-1.0000000000,-1.0000000000) 

ns 

(  0.2000001110,-0.6000000313) 

ns 

(-1.0000000000,  1.0000000000) 

T'20 

(-1 .0000000000,- 1 .0000000000) 

’21 

(  1.0000000000,-1.0000000000) 

Table  3.  Superconvergence  points  for  a  periodic  mesh  with  several  refinements:  Mas¬ 
ter  element  coordinates  of  the  points  in  the  elements  of  the  mesh  shown  in  Fig.  4c; 
bilinear  elements.  The  points  of  simultaneous  superconvergence  for  the  Xj-  and  Xj- 
derivatives  are  given.  Here  we  report  the  coordinates  of  the  superconvergence  points 
which  are  distinct  from  the  Gauss- Legendre  points.  Note  that  for  some  elements  the 
superconvergence  points  are  at  the  vertex  of  the  element. 
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Grid  with  mesh-interface  with  construned  nodes 


Bilinear  elements 

Gaoss-Legendre 

Snperconvergence  points 

Cell-size 

Points 

from  Tkble  1 

h 

wsu 

mn 

HB 

wBm 

IIB 

liB 

Element  of  size  |  at  the  mesh-interface 

0.25000 

.57517440E-I-00 

2.30 

.17657086E+00 

2.83 

0.12500 

.31086194E-I-00 

2.49 

.38090653E-01 

2.44 

.15849922E-I-00 

2.54 

.86690122E-02 

2.22 

.79755480E-01 

2.55 

.21457579E-02 

2.20 

Element  of  size 

1  one  layer  away  from  the  mesh-interface 

0.25000 

.77153390E-01 

0.31 

.77003693E-01 

1.23 

0.12500 

.34914069E-01 

0.28 

.19916182E-01 

1.28 

0.06250 

.16510365E-01 

0.27 

.52710414E-02 

1.35 

0.03125 

.77883456E-02 

0.25 

.13379941E-02 

1.37 

Table  4a.  Rate  of  convergence  at  the  superconvergence  points:  Grids  with  mesh- 
interface  with  constrained  nodes  shown  in  Figs.  6a-6d,  exact  solution  u(xi,X2)  = 

sin(7rxi)sinh(irx2),  bilinear  elements  (p  =  1).  The  values  of  at  the  Gauss- 

'  OXi  I 

Legendre  points  and  the  superconvergence  points  (from  Table  1)  are  given  in  the  squaire 

elements  of  size  —  at  the  mesh-interface  and  one  layer  away  from  the  mesh-interface 
are  given. 
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Grid  with  mesh-interface  with  constrained  nodes 

Biquadratic  elements 

Cell-size 

1 

Gauss- Legendre 

Points 

_ 

Superconvergence  points 
from  Tzble  1 

h 

1^1 

1  dll  1 

11^1 

1  0*1  1  a-L 

1^1 

1  d*i  1  mp 

1  1  l-3 

1  0*1  1  ^ 

Element  of  size  ^  at  the  mesh-interface 

0.25000 

.35616628E-01 

0.57 

1 

.34981651E-01 

2.24 

0.12500 

.45941576E-02 

0.29 

.33414467E-02 

1.71 

0.06250 

.10988098E-02 

0.28 

.41561324E-03 

1.70 

0.03125 

.26815619E-03 

0.28 

.51379346E-04 

1.69 

Element  of  size 

j  one  layer  away  from  the  mesh-interface 

0.25000 

.76959609E-02 

0.12 

.76081781E-02 

0.49 

0.12500 

.72206159E-03 

0.05 

.66641794E-03 

0.34 

0.06250 

.16100891E-03 

0.04 

.68416871E-04 

0.28 

0.03125 

.40986412E-04 

0.04 

.80145619E-05 

0.26 

Table  4b.  Rate  of  convergence  at  the  superconvergence  points:  Grids  with  mesh- 
interface  with  constrained  nodes  shown  in  Figs.  6a-6d,  exact  solution  u(X],X2)  = 

sin(7rxj)sinh(jrx2),  biquadratic  elements  (p  =  2).  The  m20cimum  values  of  at 

•  Ox  I ' 

the  Gauss- Legendre  points  and  the  superconvergence  points  (from  Table  1)  are  given 
in  the  square  elements  of  size  —  at  the  mesh-interface  and  one  layer  away  from  the 
mesh-interface  are  given. 
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Grid  with  mesh-interface  with  constrained  nodes 

Bilinear  elements 

Gauss- Legendre 

Snperconvergence  points 

Cell-size 

Points 

from  Table  1 

h 

HB 

IIB 

llB 

IIB 

Element  of  size  j  at  the  mesh-interface 

0.25000 

.95499223E-01 

1.52 

.14953978E-i-00 

2.39 

0.12500 

.20029031E-01 

1.28 

.32446341E-01 

2.08 

.44797614E-02 

1.15 

.74421844E-02 

1.91 

0.03125 

.11031333E-02 

1.13 

.18400167E-02 

1.88 

Element  of  size 

^  one  layer  away  from  the  mesh-interface 

|!iii| 

.99972051E-01 

1.60 

.99972059E-01 

1.60 

.21907125E.01 

1.40 

.20028986E-01 

1.28 

0.06250 

.47790173E-02 

1.22 

.45774021E-02 

1.17 

0.03125 

.11739024E-02 

1.20 

.11421790E-02 

1.17 

Table  5a.  Rate  of  convergence  at  the  superconvergence  points:  Grids  with  mesh- 
interface  with  constrained  nodes  shown  in  Figs.  6a-6d,  exact  solution  u{xj,X2)  = 

sin(TrXi)sinh(7ri2),  bilinear  elements  (p  =  1).  The  values  of  at  the  Gauss- 

*0X2  * 

Legendre  points  and  the  superconvergence  points  (from  Table  1)  are  given  in  the  square 


elements  of  size 
are  given. 


h 

2 


at  the  mesh-interface  and  one  layer  away  from  the  mesh-interface 
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Grid  with  mesh-interface  with  constrained  nodes 

Biquadratic  elements 

Cell-size 

Gauss- Legendre 

Points 

Superconvergence  points 
from  Table  1 

h 

m 

B 

m 

Element  of  size  ^  at  the  mesh-interface 

.29600593E-01 

1.72 

.23262232E-01 

1.49 

MnsTwriM 

.36423052E-02 

1.87 

.29870942E-02 

1.53 

0.06250 

.M509815E-03 

1.82 

.37367439E-03 

1.53 

0.03125 

.55001651E-04 

1.80 

.46708914E-04 

1.53 

Element  of  size 

j  one  layer  away  from  the  mesh-interface 

0.25000 

.52147491E-02 

0.33 

.52147491E-02 

0.33 

0.12500 

.38514170E-03 

0.20 

.38514162E-03 

0.06250 

.35773791E-04 

0.15 

.35773789E-04 

0.15 

mu 

.46018761E-05 

0.15 

.46049812E-05 

0.15 

Table  5b.  Rate  of  convergence  at  the  superconvergence  points:  Grids  with  mesh- 
interface  with  constrained  nodes  shown  in  Figs.  6a-6d,  exact  solution  u(xi,Z2)  = 

at 

CfX2 

the  Gauss- Legendre  points  and  the  superconvergence  points  (from  Table  1)  are  given 

h 

in  the  square  elements  of  size  —  at  the  mesh-interface  and  one  layer  away  from  the 
mesh-interface  are  given. 
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Grid  with  mesh-interface  with  transition  elements 

Linear  elements 

Cell-size 

Gauss-Legendre 

Points 

_ 

Superconvergence  points 
from  Table  1 

h 

1  dxi  1  g-L 

11^1  h-' 

1  dxi  1 

1^1 

1  \  „p 

m 

Element  of  size  ^  at  the  mesh-interface 

0.25000 

.26294012E-01 

0.11 

.71260207E-01 

1.14 

0.12500 

.10959010E-01 

.10239329E-01 

0.66 

0.06250 

.98340266E-02 

0.16 

.20809826E-02 

0.53 

0.03125 

.60270808E-02 

0.19 

.53221772E-03 

0.54 

Element  of  size 

^  one  layer  away  from  the  mesh-interface 

0.25000 

.24196256E-01 

0.10 

.53080555E-01 

0.85 

0.12500 

.60704882E-02 

0.05 

.74038127E-02 

0.47 

0.06250 

.16056918E-02 

0.03 

.88662264E-03 

0.23 

0.03125 

.10302041E-02 

0.03 

.23838655E-03 

0.24 

Table  6a.  Rate  of  convergence  at  the  superconvergence  points:  Grids  with  mesh- 
interface  with  transition  elements  shown  in  Figs.  7a-7d,  exact  solution  u(xi,X2)  = 

sin(7rxj)sinh(jrx2),  linear  elements  (p  =  1).  The  values  of  at  the  Gauss-Legendre 

I  OXi  I 

points  and  the  superconvergence  points  (from  Tables  2a,  2b)  are  given  in  the  square 

h 

elements  of  size  —  at  the  mesh-interface  and  one  layer  away  from  the  mesh-interface 


are  given. 
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Grid  with  mesh-interface  with  transition  elements 

Quadratic  elements 

Cell-size 

Gauss-Legendre 

Points 

Superconvergence  points 
from  Table  1 

h 

m 

B 

B 

Clement  of  size  ^  at  the  mesh-interface 

0.25000 

.40118091E-01 

0.64 

.44601274E-01 

2.86 

Mias'.  ifiTiW 

.11331178E-01 

0.73 

.11754109E-01 

6.01 

■iwii?>!w!m 

.30112286E-02 

0.77 

.10030098E-02 

4.11 

.74461581E-03 

0.77 

.12576887E-03 

4.12 

Element  of  size 

^  one  layer  away  from  the  mesh-interface 

0.25000 

.64905047E-02 

0.10 

.71044126E-02 

0.45 

0.12500 

.10139819E-02 

.11000841E-02 

0.56 

0.06250 

.24176022E-03 

.17214310E-03 

0.71 

0.03125 

.53840918E-04 

0.06 

.21691406E-04 

0.71 

Table  6b.  Rate  of  convergence  at  the  superconvergence  points:  Grids  with  mesh- 
interface  with  transition  elements  shown  in  Figs.  7a-7d,  exact  solution  u(xj,X2)  = 

sin(?rxi)sinh(?rx2))  quadratic  elements  (p  =  2).  The  maximum  values  of  -5-^  at  the 

*  UXi  * 

Gauss- Legendre  points  and  the  superconvergence  points  (from  Tables  2a,  2b)  are  given 

h 

in  the  square  elements  of  size  —  at  the  mesh-interface  and  one  layer  away  from  the 
mesh-interface  are  given. 
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Grid  with  mesh-interface  with  transition  elements 

Linear  elements 

CeU-size 

h 

Center  of  the 
element 

Superconvergence  points 
from  Table  1 

m 

IB 

Ifl 

Triangular  element  rj  at  the  mesh-interface 

.64847413E-01 

.23381462E-01 

.90781904E-02 

.43145871E-02 

.71810910E-01 

.20457913E-01 

.72098413E-02 

.18104208E-02 

Triangular  element  at  the  mesh-interface 

0.25000 

0.12500 

0.06250 

0.03125 

.55608910E-01 

.24097771E-01 

.90104544E-02 

.45067131E-02 

0.22 

0.19 

0.14 

0.14 

.64651209E-01 

.27108881E-01 

.10000891E-01 

.25104930E-02 

1.03 

1.73 

2.56 

2.57 

Table  7.  Rate  of  convergence  at  the  superconvergence  points:  Grids  with  mesh- 
interface  with  transition  elements  shown  in  Figs.  7a-7d,  exact  solution  u(xj,i2)  = 

l^-^l  are  given  at  the  center 

the  superconvergence  points 

(from  Table  2a). 


sin(?rxi)  sinh(jrx2),  linear  elements  (p  =  1).  The  values  of 
of  the  triangular  elements  T2  and  Tj  (see  Fig.  4b)  and  at 
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Rate  of  convergence  of 

Grids  irith  mesh'interface 

Bilinear  elements 

a 

A 

mm 

.38189339E-i-00 

.94191674E-01 

.24140167E-01 

6.18 

.60765580E-02 

6.22 

h 

E^h-i 

0.25000 

.42966792£-i-00 

0.12500 

.14797583E-I-00 

wmm 

.51739635E-01 

.18205665E-01 

n 

E:,h-^ 

.47927923E-I-00 

.18460718E-I-00 

.72701566E-01 

2.93 

Wt? 

.28800505E-01 

2.93 

h 

E^h-i 

IgMigJ 

.50758796E-I-00 

2.87 

.20781247E-I-00 

0.06250 

.86769758E-01 

2.78 

0.03125 

.36422194E-01 

2.78 

D 

■®io 

E,oh-^ 

0.25000 

.56625513E+00 

2.60 

0.12500 

.26052825E-I-00 

2.57 

0.06250 

.12086784E-I-00 

2.55 

0.03125 

.56320092E-01 

2.55 

h 

Eooh-' 

.62835301E-I-00 

2.51 

.31841522E-J-00 

2.51 

.15947373E-I-00 

2.55 

InRifS 

.79755480E-01 

2.55 

Table  8a.  Rate  of  convergence  of  the  discrete- norm  E^.  Meshes  shown  in  Fig.  6a- 
6d,  exact  solution  u{x^,X2)  =  sin(?rxj)sinh(7rx2),  bilinear  elements  (p  =  1).  The 
discrete-norm  was  computed  using  the  Gauss-Legendre  points  in  the  elements  in 
the  subdomain  JIq  =  (0.25,0.75)  x  (0.5,0.75)  C  fi  =  (0,1)^.  The  mesh-interface  is 
located  along  the  line  Xj  =  0.50. 


56 


Rate  of  convergence  of  Eg 

Grids  with  mesh-interface 

Biquadratic  elements 

a 

A 

Mm 

ijjjWj 

.21415864E-01 

1.37 

.27893542E-02 

1.43 

0.06250 

.35907094E-03 

1.47 

0.03125 

.45166015E-04 

1.48 

h 

.24765556E-01 

0.79 

.46104445E-02 

0.83 

.82232958E-03 

0.84 

.14526387E-03 

0.84 

h 

Ez 

E^h'i 

.27360051E-01 

.61654349E-02 

.12477526E-02 

.24915430E-03 

0.81 

a 

E4 

.29631698E-01 

.72911007E-02 

.15701595E-02 

WnSifSI 

.32847516E-03 

wbm 

L^_ 

E\o 

.35949866E-01 

.10195347E-01 

0.06250 

.25087484E-02 

0.85 

0.03125 

.58695587E-03 

0.85 

n 

Eooh-^ 

.44992746E-01 

.13862436E-01 

.34181964E-02 

R[!k“ 

.86914062E-03 

WIBM 

Table  8b.  Rate  of  convergence  of  the  discrete-norm  Eg.  Meshes  shown  in  Fig.  6a- 
6d,  exact  solution  u(ii,X2)  =  sin(7rxi)sinh(7rx2),  biquadratic  elements  (p  =  2).  The 
discrete-norm  Eg  was  computed  using  the  Gauss-Legendre  points  in  the  elements  in 
the  subdomain  Qq  =  (0.25,0.75)  x  (0.5,0.75)  C  ft  =  (0,1)^.  The  mesh-interface  is 
located  along  the  line  =  0.50. 
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Relative  error  at  the  superconvergence  points  and  Gauss  points 
in  the  elements  in  the  mesh-patch  shown  in  Fig.  8a,  p  =  1 


Element 
nomber  t 

Gauss  point 

Point  from  Table  la 

Gauss  point 

Point  from  Table  la 

Elements  at  the  mesh-interface 


10 

0.249 

11 

0.251 

12 

0.248 

13 

0.253 

14 

0.246 

15 

0.254 

16 

0.244 

0.014 

0.017 

0.015 

0.016 

0.016 

0.015 

0.017 

0.014 

0.019 

0.013 


0.084 

0.084 

0.088 

0.089 

0.089 

0.089 

0.089 

0.089 

0.088 

0.088 

0.084 

0.084 

0.079 

0.079 


Elements  one  layer  away  from  the  mesh-interface 


0.044 

0.044 

0.043 

0.043 

0.042 

0.040 

0.040 

0.039 

0.039 

0.040 

0.040 

0.042 

0.043 

0.043 

0.044 

0.044 


0.064 

0.059 

0.063 

0.058 

0.061 

0.059 

0.059 

0.059 

0.059 

0.064 

0.060 

0.066 

0.061 


0.054 

0.056 

0.044 

0.046 

0.039 


0.036 

0.036 

0.036 

0.036 

0.040 

0.039 

0.046 

0.044 

0.056 

0.054 


0.040 

0.039 

0.039 

0.040 

0.040 

0.042 

0.043 

0.043 

0.044 

0.044 
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Table  9a.  Relative  error  at  the  superconvergence  points  for  meshes  which  are  not 
locally  periodic:  Mesh  shown  in  Fig.  8a,  bilinear  elements  (p  =  1)  exact  solution 
u(£i,a:2)  =  k  =  10,  (a,  b)  =  (0,0).  The  relative  errors  in  the  Xj-  and 

Xj-derivatives  at  the  Gauss-Legendre  point  and  the  superconvergence  points  (given 
in  Table  1)  are  given  for  the  various  elements  in  Patch  I  shown  in  Fig.  8a  and  enu¬ 
merated  as  shown  in  Fig.  8b.  Note  that  the  relative  error  of  X|-derivative  at  the 
superconvergence  points  for  elements  1-16  is  around  0.02,  but  the  relative  error  at  the 
Gauss-Legendre  point  is  around  0.25. 
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Relative  error  at  the  superconvergence  points  and  Gauss  points 
in  the  elements  in  the  mesh>patch  shown  in  Fig.  8c,  p  =  1 

Element 
number  i 

Gauss  point 

Point  from  Table  la 

Gauss  point 

Point  from  Table  la 

Elements  at  the  mesh-interface 

1 

0.091 

0.084 

0.102 

2 

0.111 

0.097 

0.106 

3 

0.118 

0.101 

0.191 

4 

0.131 

0.118 

0.197 

Elements  one  layer  away  from  the  mesh-interface 

5 

0.067 

0.064 

■■ 

0.080 

6 

0.067 

0.064 

0.084 

7 

0.069 

0.064 

0.093 

8 

0.069 

0.064 

0.099 

Table  9b.  Relative  error  at  the  superconvergence  points  for  meshes  which  are  not 
locally  periodic:  Mesh  shown  in  Fig.  8c,  bilinear  elements  (p  =  1)  exact  solution 
~  ®  \  k  =  10,  (a,  b)  =  (0,0).  The  relative  errors  in  the  Xj-  and 

X2-derivatives  at  the  Gauss- Legendre  point  and  the  superconvergence  points  (given  in 
Table  1)  are  given  for  the  various  elements  in  Patch  II  shown  in  Fig.  8c  and  numbered  as 
shown  in  Fig.  8d.  Note  that  the  relative  error  of  i2-derivative  at  the  superconvergence 
points  is  smaller  than  the  relative  error  at  the  Gauss-Legendre  point. 
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Grid  with  mesh^interface  of  Fig.  4a 

Bilinear  elements 

Layer 

d 

h 

oefiB  axi 

0 

26.713 

1 

4.124 

2 

1.25 

0.454 

3 

1.75 

0.042 

Table  10a.  17%-superconvergence  in  periodic  meshes  with  local  refinements.  Grid  with 
mesh-interface  shown  in  Fig.  4a;  bilinear  elements;  Harmonic  solutions,  X|-derivative 

of  the  error.  rj%  values  at  the  (1  x  1)  Gauss-Legendre  points  in  the  elements  of  size  (— ) 

in  different  layers  away  from  the  interface.  The  distance  between  the  Gauss-Legendre 
point  considered  and  the  interface  is  denoted  by  d. 
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Grid  with  mesh'interface  of  Fig.  4a 


Biquadratic  elements 

Layer 

Point 

d 

h 

0 

1 

.1056624327 

25.137 

2 

.1056624327 

25.137 

3 

.3943375673 

1.566 

4 

.3943375673 

1.566 

1 

1 

.6056624327 

4.910 

2 

.6056624327 

4.910 

3 

.8943375673 

0.306 

4 

.8943375673 

0.306 

2 

1 

1.1056624327 

0.325 

1.1056624327 

0.325 

1.3943375673 

0.020 

1.3943375673 

0.020 

3 

1 

1.6056624327 

0.019 

2 

1.6056624327 

0.019 

3 

1.8943375673 

0.002 

4 

1.8943375673 

0.002 

Table  10b.  17%-superconvergence  in  periodic  meshes  with  local  refinements.  Grid 
with  mesh-interface  with  constrained  nodes  shown  in  Fig.  4a;  biquadratic  elements; 
Harmonic  solutions,  Xj-derivative  of  the  error.  r}%  values  at  the  (2  x  2)  Gauss- Legendre 

points  in  the  elements  of  size  (— )  in  different  layers  away  from  the  interface.  (The 

master  element  coordinates  of  the  points  in  the  element  are:  Point  1  =  {—p,  ~7=); 

v3  v3 

Point  2  =  (^,  ^);  Point  3  =  (^,  ^);  Point  4  =  -^)).  The  distance 

between  the  Gauss-Legendre  point  considered  and  the  interf£M:e  is  denoted  by  d. 


Table  10c.  t7%-supercoiivergence  in  periodic  meshes  with  local  refinements.  Grid  with 
mesh-interface  with  transition  elements  shown  in  Fig.  4b;  linear  elements;  Harmonic 


solutions,  X] -derivative  of  the  error.  ti%  values  at  the  (1x1)  Gauss- Legendre  points 

in  the  elements  of  size  (2)  different  layers  away  from  the  interface.  The  distance 
between  the  Gauss-Legendre  point  considered  and  the  interface  is  denoted  by  d. 


Grid  with  mesh-interface  of  Fig.  4b 

Quadratic  elements 

Layer 

Point 

d 

h 

a 

_  max  Q, w,  1, f) 

0X1 

0 

1 

.1056624327 

36.202 

2 

.1056624327 

16.369 

3 

.3043375673 

6.812 

4 

.3943375673 

7.243 

1 

1 

.6056624327 

25.589 

.6056624327 

12.789 

.8943375673 

0.813 

.8943375673 

1.544 

2 

1 

1.1056624327 

15.290 

2 

1.1056624327 

5.707 

3 

1.3943375673 

0.546 

4 

1.3943375673 

0.003 

3 

1 

1.6056624327 

3.001 

2 

1.6056624327 

0.004 

3 

1.8943375673 

0.000 

4 

1.8943375673 

0.000 

Table  lOd.  77%-superconvergence  in  periodic  meshes  with  local  refinements.  Grid 
with  mesh-interface  with  transition  elements  shown  in  Fig.  4b;  quadratic  elements; 
Harmonic  solutions,  Xj-derivative  of  the  error.  values  at  the  (2  x  2)  Gauss- Legendre 

points  in  the  elements  of  size  (— )  in  different  layers  away  from  the  interface.  (The 

-1  -1 

master  element  coordinates  of  the  points  in  the  element  are:  Point  1  = 

Point  2  =  ^);  Point  3  =  Point  4  =  (^,  ^)).  The  distance 

between  the  Gauss-Legendre  point  considered  and  the  interface  is  denoted  by  d. 
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